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ABSTRACT 

We present a sample of 290 24/im-selected active galactic nuclei (AGNs) mostly 
at z ~ 0.3 — 2.5, within 5.2 deg 2 distributed as 25' x 25' fields around each of 30 
galaxy clusters in the Local Cluster Substructure Survey (LoCuSS). The sample is nearly 
complete to 1 mJy at 24/jm, and has a rich multi-wavelength set of ancillary data; 162 
are detected by Herschel. We use spectral templates for AGNs, stellar populations, and 
infrared emission by star forming galaxies to decompose the spectral energy distributions 
(SEDs) of these AGNs and their host galaxies, and estimate their star formation rates 
(SFRs), AGN luminosities, and host galaxy stellar masses. The set of templates is 
relatively simple: a standard Type-1 quasar template; another for the photospheric 
output of the stellar population; and a far infrared star-formng template. For the 
Type-2 AGN SEDs, we substitute templates including internal obscuration, and some 
Type-1 objects require a warm component (T > 50 K). The individually Herschel- 
detected Type-1 AGNs and a subset of 17 Type-2 ones typically have luminosities 
> 10 45 ergs s' 1 , and supermassive black holes of ~ 3 x 10 8 Mq emitting at ~ 10% of 
the Eddington rate. We find them in about twice the numbers of AGN identified in 
SDSS data in the same fields, i.e., they represent typical high luminosity AGN, not an 
infrared-selected minority. These AGNs and their host galaxies are studied further in 
an accompanying paper. 

Subject headings: galaxies: active-quasars: general-infrared: galaxies 


1. Introduction 

The bright continua of active galactic nuclei (AGNs) in the X-ray, UV and optical are powered 
directly by accretion - the growth of super massive black holes (SMBHs). At the current epoch, 
there is a tight correlation between the SMBH masses and the host galaxy stellar bulge masses, 
indicating a link between the integrated accretion by black holes and the star formation in their host 
galaxies (e.g., Magorrian et al. 1998; Tremaine et al. 2002). The level of accretion is indicated by 
a variety of metrics, e.g. X-rays, optical emission lines, and optical-IR continua. These indicators 
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can drown out many metrics for the level of star formation, but it is thought that the far infrared 
(FIR) emission remains dominated by this mechanism, providing a strong motivation for studies of 
the FIR outputs of galaxies with active nuclei. 


Prior to Herschel, the measurements of rest-frame FIR emission from luminous AGNs were 
limited to a small population (e.g., Omont et al. 2001; Haas et al. 2003; Dicken et al. 2008). With 
the advent of Herschel , it is possible to study the FIR properties efficiently for a large sample of 


AGNs (e.g., 

Hatziminaoglou et al. 

2010 

Shao et al. 

2010 

Mullaney et al. 

2012 

Rosario et al. 

2013, Leipski et al. 

2013, 

2014 

). To augment these studies, we describe Herschel measurements of 


205 Type 1 AGNs uniformly selected from a 5.2 deg 2 survey area; the sample is nearly complete at 
24 fim down to 1 mJy, and verified by spectroscopy. These AGNs are complemented by 85 Type 
2 objects similarly selected from a 3.6 deg 2 subset of the same data. A multi-band dataset from 
the UV to the FIR, including optical spectroscopy, allows detailed study of the AGNs and their 
host galaxies in this sample. We use spectral templates for 1.) the UV to far infrared output of 
AGNs; 2.) stellar populations; and 3.) the infrared emission of star forming galaxies to decompose 
their spectral energy distributions (SEDs) and estimate their IR star formation luminosities, AGN 
luminosities, and their host galaxy stellar masses. The black hole masses are estimated from the 
Type-1 AGN broad optical emission lines, and from the stellar masses for the Type-2 objects. More 
than 55% of the sample are detected individually by Herschel, and we stack the signals from the 
rest for comparison purposes. This paper presents the data and a basic analysis; in a companion 
paper (Xu et al. 2015), we use these results to explore the evolutionary stage of the AGN-host 
galaxies - whether they are starbursting or normal star-forming galaxies; and whether there is a 
causal connection between nuclear and SF activity for these objects. 

This paper is structured as follows. In Section 2 we present the data, and in Section 3 we de¬ 
scribe the selection and completeness of our sample of Type-1 AGNs. We also show their SEDs and 
discuss the FIR excesses. Section 4 is a parallel discussion of the Type-2 objects. In Section 5 we 
analyze both samples together, calculating the physical properties of the AGNs and their host galax¬ 
ies, such as Eddington rates, host stellar masses and star formation rates. A summary is provided 
in Section 6. Throughout this paper we assume Dm = 0.3, Da = 0.7, and Hq = 70kms _1 Mpc^ 1 . 

2. DATA 

2.1. LIRAS: LoCuSS Infrared AGN Survey 

The Local Cluster Substructure Survey (LoCuSS^ is a large survey of X-ray luminous galaxy 
clusters at z = 0.15 — 0.3 (e.g., Smith et al. 2010). This paper exploits the extensive LoCuSS 
multi-wavelength data set for 30 clusters, which includes data from Chandra, GALEX, SUBARU, 
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UKIRT, Spitzer/ MIPS, and Herschel. The Spitzer and Herschel data cover a total area of ~ 5.2 
deg 2 (25 / x 25 ; x 30), at the central coordinates listed in Table [l] Since most cluster members are 
members of the old galaxy population, which is not bright in the MIR and FIR, this wide-held 
coverage allows us to conduct a serendipitous AGN survey independent of the existence of galaxy 
clusters in the observed fields. In LIRAS (LocuSS Infrared AGN Survey), we take advantage of 
these multi-wavelength datasets (described in this section) to study the properties of a 24 ^in- 
selected IR luminous Type-1 AGN sample over the entire 5.2 deg 2 area (as discussed in Section 3 
below). In Section 4, we show how we also identified Type-2 AGNs selected from 21 out of the 30 
cluster fields (i.e., 3.6 deg 2 ), as indicated in Table [lj 


2.2. Mid-infrared Observations 


Each cluster held was observed at 24 ^m between November 2007 and November 2008 with 
MIPS ( Rieke, et al.|2004 ) on Spitzer ( Werner et al.|2004 ), utilizing a 5 X 5 grid of pointings in hxed 
cluster or raster mode (PID: 40872; PI: G.P. Smith). Two cycles of small-held photometry with a 
frame time of 3 s were performed at each grid point, for a total per pixel exposure time of 90 s. 
The central 5' x 5' of some clusters had already been imaged by GTO program 83 to greater depth 
(~ 3000 s/pixel). All the available data were combined for our survey. The images were processed 
with the MIPS Data Analysis Tool (DAT; Gordon et al. 2005). The beam size at 24 pm is 5 / . / 9 with 
2 / . , 49 pixels; the images were combined with a pixel scale of T/245, half the physical pixel scale. The 
24 pm fluxes were measured by SExtractor (Bertin Sz Arno uts|l9 96) within a hxed circular aperture 
of diameter 21” and with an aperture correction of a factor 1.29. The 90% completeness limits at 
24 pm are in the range 300-500 pJy. Details of the reduction, source extraction and photometry 
can be found in ?. 


WISI0 data are also available in our survey helds. We utilize WISE 3.4 pm, 4.6 pm, and 12 
pm measurements in our SED decomposition htting. The detection limit at 22 pm is 6 mJy, an 
order of magnitude higher than achieved with MIPS. Since this band is so close to the MIPS 24 
pm one, we do not use it. 


2.3. Far-infrared Data 

Our Herschel (Pilbratt et al. 2010) data were taken between 22 December 2009 and 10 October 
2011 (LoCuSS Herschel key Programme, Smith et al. 2010). Each cluster held was observed with 
both the Photodetector Array Camera and Spectrometer (PACS; Poglitsch et al. 2010) at 100 
and 160 pm, and the Spectral and Photometric Imaging Receiver (SPIRE; Griffin et al. 2010) 
at 250, 350, and 500 pm. The images were reduced using HIPE V6.0 (Ott 2010). Because the 


2 http://irsa.ipac.caltech.edu/Missions/wise.html 






-4- 


PACS images are relatively shallow and the beam size is relatively small, confusion is not an issue 
and the photometry could be performed with SExtractor. However, confusion is an issue for the 
SPIRE data. Therefore the photometry of the SPIRE images was performed with IRAF/DAOphot, 
using the 24 fi m source positions (for all sources above the 3er detection limit) as priors to position 
the PSF (point spread function) on the SPIRE maps. We rotated the Herschel PSF to match the 
position angle of each map, registered the 24 /xrn and Herschel maps with the isolated point sources, 
and then fixed the source positions. We extracted fluxes on the SPIRE maps using the empirical 
fine scale PSF provided by the HSC[^] instead of constructing one from our own data, because of the 
lack of isolated point sources with high ratios of signal to noise on our maps. Parameters adopted 
for the maps and photometry (such as the pixel size, FWHM of point source, photometry aperture 
radius, aperture correction, and sky annulus radius) are summarized in Table [2j 


2.4. Near-infrared, Optical, and Ultraviolet Data 


Near-infrared images of 26 of the 30 cluster fields were obtained with WFCAM (Casali et al. 
2007) at J- and K-bands on the 3.8-m United Kingdom Infrared Telescope (UKIRT) in service 
mode over multiple semesters starting in March 2008. The data acquisition used the same strategy 
as was used by the UKIDSS Deep Extragalactic Survey (Lawrence et al. 2007), covering 52' x 52' to 
depths of J ~ 21, K ~ 19, with exposure times of 640 s, pixel size of 0 ".2 (half the physical pixel size) 
and PSF FWHMs ~ if.7 —1".‘2. The remaining four cluster fields were observed with NEWFIRM on 
the 4.0-m Mayall telescope at Kitt Peak on 17 May 2008 and 28 December 2008. The NEWFIRM 
data consist of dithered and stacked J- and K-band images covering fields of 27 ; x 27 / with a 0 ".4 
pixel-scale and PSF FWHM ~ l'.'O — T/5. The total exposure times in each filter were 1800 s, and 
the images also reach depths of J~21 mag and K~19 mag (Vega, 5-sigma). 


SDSS photometry (Data Release 7) is available for 26 out of 30 cluster fields, covering a 
total of 4.51 deg 2 survey area. The SDSS five band photometry we used is corrected for Galactic 


extinction. Optical images in R or I band using Subaru/Suprime-Cam (Okabe et al. 2010) with 
seeing ~ 0 / . , 6 allow us to study the morphology of the Type 2 AGN hosts. The data were reduced as 
described by Okabe &: Umetsu (2008), using the Suprime-Cam pipeline software SDFRED for flat- 
fielding, instrumental distortion correction, differential refraction, PSF matching, sky subtraction, 
and stacking. The astrometric solution was based on 2MASS stars. Standard stars were interspersed 
with the cluster imaging. Further information about the optical imaging, including the initial 


publication of the data for most of the clusters, can be found in Okabe et al. (2010) 


GALEX NUV observations were obtained for 26 of the cluster fields (omitting those for A586, 
A689, A2485, and RXJ0142), and simultaneously in the FUV for 21 fields under Guest Investigator 
Programs GI4-090 and GI6-046 (Pis G. P. Smith and S. Moran, respectively). The exposure times 
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ranged from 3 to 29 ks; additional details about these observations and their reduction can be 
found in Haines et al. (2015). 


2.5. Chandra X-ray imaging 


Twenty one of the 30 clusters were observed with Chandra in the I mode of the Advanced 
Camera for Imaging Spectroscopy (ACIS-I), which has a field of view of 16(9 x 16(9. Seven more 
were observed with the ACIS-S (8(3 x 8(3 FOV). Two of the clusters (Abell 2345 and Abell 291) do 
not have X-ray data. The exposure times for the cluster fields range from 10 ks to 100 ks (Table 1 
in Haines et al. 2012), with a typical integration time of 20 ks. 


Sanderson et al. (2009) discuss the reduction and analysis of the X-ray observations. To 
detect X-ray point-sources that are potential AGNs, we used the wavelet-detection algorithm CIAO 
WAVDETECT; a minimum of six counts in the broad energy (0.3-7 keV) range was the threshold for 
source detection. The observations in this band were converted to fluxes assuming a T=1.7 power- 
law spectrum with Galactic absorption, following Renter et al. (2005). The X-ray flux sensitivity 
limit for the cluster fields ranges from 6 x 10~ 16 to 8 X 10~ 15 ergs cnr 2 s -1 with a median of 
3.5 x 10 -1 ' 5 ergs cnr’ s -1 . We calculated the luminosities for all sources with redshifts, assuming 


,—2 0 -l 

VVt L 

K-corrections of the form (1 + z) l ~ 2 . 
detected AGN. 


We also calculated the X-ray luminosity limit for each non- 


2.6. Spectroscopic data 


We use spectra from ACReS (the Arizona Cluster Redshift Survey; Pereira et al. 2015 in 
preparation) a long-term spectroscopic program to observe the fields of the 30 galaxy clusters with 
MMT/Hectospec. Hectospec is a 300-fiber multi-object spectrograph with a circular field of view of 
1° diameter (Fabricant et al. 2005) and fibers that project to T/5 on the sky, mounted on the 6.5-rn 
MMT at Mount Hopkins, Arizona. We used the 270 line grating, which provides a wide wavelength 
range (3650A-9200A) at 6.2 A resolution (R = 1000). This ensures coverage of the most important 
emission lines suitable for identifying AGNs. The spectroscopy data were reduced using HSREE(^J 
Redshifts were determined by comparison of the reduced spectra with stellar, galaxy and quasar 
template spectra, choosing the template and redshift that minimized the y 2 between model and 


data. The target selection is described in detail in Haines et al. (2013). Virtually all sources with 


24 /im flux above 1 mJy, and K-band < 19 mag, were targeted by Hectospec. See Section [3] and 
Appendix [B] for a summary of the spectroscopic coverage of the 24 /rrn sources. 


4 http://mmto.org/ rcool/ 
















3. 24 /tm-Selected Type-1 AGNs 


The mid-IR continuum emission of Type-1 AGNs arises from warm dust heated by the AGN 
(e.g., Rieke 1978; Polletta et al. 2000; Haas et al. 2003), and on average there are only modest 
variations among quasars in the average fraction of the bolometric luminosity emitted at these 
wavelengths (Krawczyk et al. 2013). Therefore, the 24 /ini selection of Type-1 AGN is expected 
to be highly efficient and complete. 


3.1. Type-1 AGN Identification 

3.1.1. Approach 

Members of our sample of LIRAS Type-1 AGNs (see Table 3 for those detected and Table 4 
for those undetected with Herschel) were required to have: 

1. Spitzer /MIPS 24 /mi flux densities above 1 mJy; and 

2. optical spectra showing broad emission lines with full width at half maximum (FWHM) 
> 1200 km s -1 . 


All sources in the surveyed regions with 24 /mi flux densities above 1 mJy and with K-band 
< 19 mag were given the highest priority for spectroscopy, irrespective of near-IR color or morphol¬ 
ogy (resolved/unresolved in the near-IR). We excluded those that were clearly stars (being both 
unresolved in the K-band data, and having blue near-IR colors (J — K < 1.0)) and asteroids (very 
luminous at 24 /mi, but with no counterparts in all other bands). Over the 5.2 deg 2 survey area 
covered at both 24 /mi and in the Herschel bands, there were 2439 sources with 24 /mi flux density 
above 1 mJy, of which 1827 remained after excluding stars, asteroids, and sources with no opti- 
cal/NIR counterparts. From this list, 1729 sources were observed with Hectospec while another 18 
have spectra from SDSS. Therefore, the completeness of spectroscopic coverage is about 94.6% (See 
Appendix A for a summary of the reasons that we missed 5.4% of the targets). Among the 1729 
sources targeted by Hectospec, 1263 yielded spectroscopic redshifts with a corresponding success 
rate of 73%. 


To identify Type-1 AGNs, we fitted each emission line in the optical spectra with single or 
double Gaussian profiles. We list the FWHM of typical broad emission lines in Table 5. Sources 
showing emission lines (specifically, Mgn, Civ, H/3, or Hcc) with FWHM over 1200 km s -1 were 
selected as Type-1 AGNs (H ao et ah] 2005). Finally 205 sources satisfied our Type-1 AGN selection 
criteria, 177 confirmed with Hectospec and 28 confirmed with SDSS^j More details of the sample 


s The spectra were inspected visually to confirm the classification. One ambiguous case (#6 J084352.28+292854.0) 
was retained because its SED decomposition fit (see Section 3.4) supported the classification. J084234.94+362503.2 
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selection can be found in Appendix A. Figure [l] shows the 24 /xrn flux distribution and spectroscopic 
status of the members of our sample. Figure [2] shows the redshift distribution of the AGNs; virtually 
all of them are far behind the clusters. They are also far from the cluster centers (typically by 1' or 
more), so they are not significantly magnified. We exclude any AGNs in the clusters according to 
the source redshift. 


3.1.2. Results 

Of the 205 24 ^urn-selected Type-1 AGNs we identified, 107 are securely detected by Herschel, 
101 in at least two bands. For these sources, Table [3] presents coordinates, redshifts, and observed 
flux densities in the near-infrared (J and I\ bands from UKIRT and NEWFIRM), mid-infrared 
(3.4, 4.6, and 12 ^xm from WISE; 24 fim from Spitzer), and far-infrared (100 /xm, 160 /xm, 250 
H m, 350 n m, and 500 ^m from Herschel). Similar information is provided in Table 4 for the 
sources not detected with Herschel. Below we discuss the basic properties of the Herschel- detected 
sources and compare them with those not detected by Herschel through a stacking analysis of the 
latter. In Section 3.5, we show that the SDSS colors of the Herschel- detected and Herschel- non- 
detected sample members do not differ significantly, showing that the dust emitting in the FIR 
is not producing significant extinction along the line of sight toward the AGNs. This result is 
consistent with the hypothesis that this dust is not associated with the active nuclei and simplifies 
the comparison of the Herschel detected and non-detected sources. 


3.1.3. Completeness 


We ran simple simulations to test the completeness of the Type-1 AGN identification. The 
two plots in the upper panel of Figure [3] show the expected apparent K-band and r' band (SDSS) 
magnitudes as a function of 24 p.m flux density for Type-1 AGNs, represented by the template of 


Elvis et al. (1994). The K and r' band magnitudes are functions of redshift and dust extinction. To 


simulate these effects, the Elvis et al. (1994) template was redshifted incrementally over the range 
z = 0 to z = 3.6, covering the redshift range over which we identified AGNs. We added extinction 


using a composite reddening law: 1.) the Galactic extinction curve above 1 /im (Rieke &: Lebofsky 


1985 

); and 2.) a SMC extinction curve below 1 /xm ( 

Gordon et al. 

2003 

6 


the value of Ay ranged 

from 0 to 1.5. We rescaled the AGN template to make the 24 fi m flux density always above 1 mJy. 
We used the same range of the flux level and Ay to generate the data points in the two plots in the 


(#13) falls slightly below our line width criterion but was retained because spectropolarimetry shows it to have a 
hidden Type 1 nucleus (Zakamska et al. 2005). 

6 Studies show that the reddening toward quasars is dominated by SMC-like dust at the quasar redshift (e.g. 
Richards et al. 2003; Hopkins, et al. 2004; also see the “Gray” extinction curve in Czerny et al. 2004 and Gaskell et 
al. 2004) 













upper panel of Figure [3j in other words, the data points in these two plots are from the same AGN 
populations. The simulation then lets us compare the incompleteness resulting from our detection 
limits for K-band and for optical spectroscopy. 

Sources fainter than our K-band detection limit of 19 mag (Vega, 5-sigma) were not targeted 
in ACReS, but this K-band limit does not affect our AGN survey significantly. As shown in the 
upper left of Figure [3j all sources with 24 /xrn flux density above 1 mJy and Ay smaller than 1.5 
mag have K-band magnitudes brighter than 19 (Vega), and therefore were targeted. The K-band 
limit would only exclude a few very red AGNs close to the 24 gm 1 mJy flux density cutoff. The r' 
band simulation shows that such very red AGNs would drop below the limit for successful optical 
spectra. As can be seen from the lower left of Figure [3j the K-band distribution of our 205 Type- 
1 AGNs declines steeply from 17.5 to 19 mag (Vega), probably because the optical spectroscopy 
sets a tighter constraint for inclusion in our sample than does the K-band cutoff. As expected, 
reddening affects the r' band more than the K-band (see upper panels of Figure [3]). The lower right 
of Figure |3] compares the distribution of r' magnitude for targets that were put on Hectospec fibers 
and targets where emission lines were detected (for sources with 24 gm flux density above 1 mJy). 
The redshift success rate starts declining from r' = 19.5 mag (AB), and declines steeply for sources 
fainter than 20.5 mag (AB). We conclude that our spectroscopic survey is incomplete at the lowest 
24 /in 1 flux density levels for red sources with Ay above 0.5 mag. 

However, Figure [3] implies that the completeness will increase rapidly as the 24 gm threshold 
is raised above 1 mJy. In confirmation, Figure [l] shows that the fraction of Type-1 AGNs within 
the total sample targeted for spectra is roughly constant above 2 mJy, but is somewhat lower in 
the 1—2 mJy bin. This drop toward 1 mJy is the behavior expected from incompleteness, but a 
part of the drop is also likely to be intrinsic, as shown by Brand et al. (2006). From Brand et al. 
(Figure 4, both data and Pearson models), we estimate that the intrinsic fraction at 1 — 1.5 mJy 
should be about 70% of the asymptotic value at larger flux densities. We then predict 100 AGNs 
in the 1 — 1.5 mJy bin, where only 68 are detected, i.e., we are potentially missing about 30 AGNs. 
There is no evidence from the counts for any missing AGNs in the 1.5 — 2 mJy or higher bins. We 
therefore estimate that the incompleteness in our sample is about 30/235 = 13%, concentrated in 
the 1 — 1.5 mJy range and largely due to incompleteness in the optical spectroscopy. Combining 
with the incompleteness of 5.4% in the spectroscopy itself, the total of missing sources is about 
18%. 


In addition to the missing Type-1 AGNs, our sample will miss other types of active object. 
For example, from Dey et al. (2008), nearly half of the objects not targeted for spectroscopy 
because they did not have optical counterparts (see Appendix A) are likely to be dust-obscured 
galaxies. Because of their low accretion efficiency and low Eddington ratios, our sample will also not 
include a significant number of jet-mode AGN (Yuan & Narayan 2014). Our sample is therefore 
confined to traditional Type-1 AGN identified by optical spectroscopy, selected to a uniform level 
of mid-infrared flux density. 
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3.2. Type 1 AGN Properties 

3.2.1. Redshift Distribution 

Figure [ 2 ] shows the redshift distribution of our Type-1 AGN samplcJ^J We omit four sources 
at z > 3 (numbers 104 - 107 in Table 3) from further analysis because of the very small number 
statistics in our sample at these redshifts. The distributions for the Herschel- detected subsample 
and the Herschel- non-detected one are very similar. The two-sanrple Kolmogorov-Smirnov test 
(K-S test) is consistent at the 5% level (P-value= 0.053) with the null hypothesis that these two 
subsamples are drawn from the same distribution; that is, they are statistically indistinguishable. 
This situation is only possible if the luminosity of the FIR excess (which we will attribute to star 
formation) grows with the increasing AGN luminosity that results from the flux limit of our AGN 
selection. 


3.2.2. Spectral Energy Distributions 

Three typical examples of the SEDs of these AGNs are illustrated in Figure[4j In vf v — A units, 
some of the SEDs look flat from the optical to the FIR, while some show peaks in the UV and 
optical, or in the FIR, or both. The FIR peak has a Rayleigh-Jeans tail, declining steeply toward 
the mm-wave. The SEDs of some sources at lower redshift show a peak near 1 /nn (rest-frame). 

By stacking signals over a large number of source positions, we can study the far-infrared 
properties of sources too faint to detect with Herschel individually. We stacked the signals for 
Herschel- non-detected Type-1 AGNs in three redshift bins: 0.1-0.7, 0.7-1.2, and 1.2-1.9. There 
are 24, 24, and 18 sources in these three bins, respectively. All the sources were well-detected from 
the UV to 24 ^m. At these wavelengths, we took the average of the measured fluxes in each band 
for all of the sources in a redshift bin, after eliminating the 3-sigma outliers. For the five Herschel 
bands, we registered PACS/SPIRE images with the 24 /xm image by aligning the bright isolated 
point sources, and then isolated a small image centered on the source position. We checked the 
images individually, and rejected those contaminated by close bright objects. For each redshift bin 
and each Herschel band, we then clipped the 3-sigma outliers for each aligned pixel over the full set 
of images. Since the Herschel maps were roughly uniform in exposure and noise, we could then take 
the straight average of all the remaining data. Where there was a detection on the stacked image, 
we did aperture photometry, including applying aperture corrections according to the parameters 
in Table [ 2 ] The resulting average SEDs in the three redshift bins (i.e., z = 0.1-0.7, 0.7-1.2, and 
1.2-1.9) are shown in Figure [5j to be compared with the SEDs of the sources detected individually 
(e.g., Figure |4j). The SEDs are similar, except the FIR peak is modestly weaker in the SEDs of the 


7 From SDSS spectroscopy, there are a number of sources detected by Herschel at z > 2; however, they are not 
included in our study because at these redshifts the Herschel bands may be significantly influenced by emission from 
the AGN. 
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stacked signals. 


3.2.3. Nature of the Far-Infrared Excess 


We now explore the nature of the far infrared peak in the quasar SEDs. We display the [250/24 
fim] flux ratio for the individual Herschel -detected objects as a function of redshift in Figure [6j 
along with the upper limits for the Herschel- non-detected AGNs and the three average values for 
the stacked SEDs. We also plot the flux ratios of a theoretical AGN template from Fritz et al. 


(2006) and a Type-1 AGN template from Elvis et al. (1994). The individual sources, upper limits, 


and stacked points all fall far above the model predictions. That is, AGN dust heating is unlikely 


to produce adequate FIR emission even if large (kpc-scale) tori are assumed (Fritz et al. 2006 


Ballantyne et al. 2006). Therefore, the prominent far-infrared excess over the AGN templates 


indicates that a star formation component may contribute significantly to the FIR. Support for 
this conclusion is provided in Section 3.4, where we show that the far infrared spectral energy 
distributions for most of the sources are similar to those of normal star-forming galaxies of similar 
luminosity. 


3.3. The FIR Dust Temperature and Mass 


We estimate the temperature of the FIR-emitting dust using a single-temperature grey-body 
model, of the form B U (T^ j)[l — e ~ Td ], where B U (T) is the Planck function and Td is the frequency- 
dependent dust optical depth. The dust is optically thin in the FIR, and we have 1 — e~ Td « r{y) = 
t{vq)(v / vq)P . Studies of local galaxies (Hildebrand |l983 Dunne fc Eales||2001 Gordon et al. 


show that a value of j3 = 1.5 is a good estimate of the emissivity index for active star formation 


regions. We use the same criteria as in Hwang et al. (2010) to select sources with well-sampled SEDs 


around the peak of the FIR emission, i.e., there should be at least one flux measurement shortwards 
and longwards of the FIR peak, and the FIR SED should be physical (convex, not concave). There 
are 36 Type-1 AGNs in our sample that meet these conditions. The dust temperatures of these 
AGNs have large scatter from 22 K to 62 K. As shown in Figure [7j we compared the results with 
the luminosity-temperature relation for star-forming galaxies at z = 0.1 — 2 derived from Hei'MES 


and PEP data in the COSMOS, GOODS-S and GOODS-N fields (Symeonidis et al. 2013). The 


majority of the AGNs with cold dust temperatures lie within the 1-sigma range of dust temperature 
of star-forming galaxies, suggesting that the origin of their FIR emission is the same as for the cold 
dust in normal galaxies, that is star formation. Most such AGNs are at the lower luminosity end. 
The galaxies with temperatures significantly above expectations for star-forming galaxies all have 


warm components (Section 3.4.3) that make determining the behavior of any cold dust ambiguous. 


For the galaxies where the FIR is dominated by the emission of cold dust, we can estimate the 
required mass of interstellar material. In the optically thin limit, 
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(1) 


where S u is the observed flux density at v , u r is the rest-frame frequency, and the mass absorption 
coefficient, K(y r ) = ko(d/D))^, is approximated by a power law. Here we take fto(125 //m) = 
18.75 cm 2 g _1 from Hildebrand (1983). The FIR-emitting dust mass ranges up to 9 x 10 8 M 0 , 


with a median value of 2 x 10 s M 0 , indicating huge reservoirs of gas in these systems. The large 
FIR luminosities for most of the other Herschel -detected systems also indicate large amounts of 
interstellar material. 


3.4. Type-1 SED Decomposition 

We will show that the SEDs in Figures [4] and [5] can be explained in terms of three dominant 
SED components; 1.) an averaged Type-1 quasar continuum to supply the UV (the “big blue 
bump” and to fill in the near-infrared; 2.) the far infrared SED of a luminous star-forming galaxy; 
and 3.) the SED of a moderately old stellar population, which peaks at wavelengths slightly longer 
than 1 /irn. These components are shown in the figures. In a few cases we need to add a warm far 
infrared component. The overall simplicity of the SED (only three dominant components) underlies 
the success of our fitting it to extract the underlying properties of the sources. SED decomposition 
can determine quantitatively the relative contribution of the AGN and the old stellar population 
in the NIR and of the AGN and star formation in the FIR. 


3-4-1. Decomposition Procedures 


One of the uncertainties for the SED decomposition is the lack of a clean template of a naked 
Type-1 AGN SED. Numerical AGN models and semi-analytic models provide candidate templates. 
The numerical models assume a central point-like energy source with a broken power-law SED 
surrounded by a smooth or clumpy dust distribution, and then solve the radiative transfer equation 
(e.g. Fritz et al. 2006). Templates generated using this method must make assumptions about 
the dust distribution geometry and compositions. For semi-analytic methods (e.g. Mullaney et al. 
2011; Sajina et al. 2012), the SED is taken to be a broken power-law based on physical assumptions 
for hot and warm components, and a modified blackbody beyond a given wavelength. Thus, both 
the numerical and semi-analytic methods are based on assumptions and introduce a number of free 
parameters. 


To minimize the number of free parameters in our fits, we use an empirical AGN SED template 
(Elvis et al. 1994) to determine the far-IR properties of our sources, and to estimate the relative 
contribution from the AGN and the host. This template may include a contribution from star 
formation and hence be too bright in the far infrared. In Appendix C, we derive a correction to 
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the template that represents a bounding condition for maximal star formation. We begin with the 


published template and consider the implications of this second case in Section 3.4.2 


We model the observed SED as a linear sum of a stellar component, a star formation com¬ 
ponent, and an AGN component. We also assume the UV emission arises from the AGN rather 
than star formation. The star formation and AGN are taken to be independent and not to affect 


each other. For the star formation component, we use the 10 infrared galaxy templates from Rieke 
et al. (2009) for luminosities of 10 9 ' 75 Lq to 10 12 Lq. This restricted range of IR luminosity has 


been shown to give appropriate FIR SEDs for galaxies at the redshifts in our study (Rujopakarn 


et al. 2013). For the stellar component, we use 24 simple stellar population SEDs from Bruzual 


and Chariot (2003), assuming a Salpeter IMF, Padova evolutionary tracks, and solar metallicity, 
at ages from 0.4 Myrs to 13 Gyrs. Populations older than 1 Gyr have very similar SEDs in the 
NIR and beyond, with any differences confined to shorter wavelengths (where the stellar output is 
overwhelmed by that of the AGN). We also make sure the stellar component is not older than the 
age of the Universe. As necessary, we apply foreground extinction to the AGN template; we use the 
Galactic extinction curve above 1 /mi (Rieke & Lebofsky 1985), and use a SMC extinction curve 
below 1 /urn (Gordon et al. 2003) (See Section 3.1.3 for more details about the reddening model). 


We fit the above model to the data for rest wavelengths from 1216 A to 1000 /irn (we do 
not use the photometry short of 1216 A due to Lyman a forest absorption). Specifically, we 
use photometry from GALEX, SDSS, J, K, WISE (3.4 /mi, 4.6 /mi, and 12 /mi bands), 24 /mi, 
Herschel (100 /mi, 160 /mi, 250 /mi, 350 /mi, and 500 /mi bands) if they are available. There are 
six degrees of freedom in the fitting: choosing the best fitting (1) infrared star formation and (2) 
optical/near infrared stellar population templates from the libraries, (3) determining the extinction 
to the AGN, and normalizing (4) the Type-1 AGN, (5) star formation, and (6) stellar population 
templates. We use Levenberg-Marquardt least-squares fitting to find the best solution among these 
degrees of freedom. 


3-4-2. SED Decomposition Results 

Figures [4] and [5] display examples of the SED decomposition results. In the rest-frame UV 
band and the MIR, the AGN always dominates the emission for our sample. The UV and optical 
are dominated by thermal emission from an accretion disk, and a significant portion of the NIR 
output is from hot dust warmed by the AGN. Many SEDs show a minimum at 1 /mi (rest-frame), 
resulting from the upper temperature limit (sublimation temperature) for grains that can survive in 
the vicinity of an AGN. The stellar component contributes in the NIR, sometimes producing a NIR 
peak, or making the SED flat near 1 /mi. The contribution from the stellar component is generally 
more significant at lower redshift. It fades quickly as the redshifts of the sources increase, as a result 
of our 24 /mi selection being dominated by the AGN and selecting increasingly luminous AGNs 
with increasing redshift. A star formation component is needed for 95% (102 out of 107) of the 
Herschel -detected sources, but the contribution of this component in the FIR varies substantially 
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from source to source. A star formation component is also required in all three redshift bins for 
the stacked SEDs of the sources not detected individually with Herschel. 

We have tested the necessity for star formation in the SED fits on the minority of sources where 
the FIR luminosity is relatively small, namely quasars number 8, 24, 27, 38, 48, 49, 61, 64, 69, 71, 
77, 90 (distributed over nearly the full redshift range of the sample, i.e., from z = 0.4 to z = 2.1). 
We determined the minimum y 2 for fits to the measurements at rest-frame wavelengths > 6 fi m, 
assuming 20% minimum effective error (larger errors for low signal to noise measurements) for each 
photometric measurement and using just the Elvis template. We evaluated the quality of the fits 
based on the values of x 2 and the number of degrees of freedom for each galaxy, and then compared 
with the result with a star forming template added. The probability of the fit being adequate with 
the Elvis template alone was < 0.3% in every case. With the star forming template added, the 
probability that the fit was adequate was > 15% for seven (of twelve) cases, > 1% in three more, 
and was always much larger than the probability without the star forming template. The two cases 
with bad fits had far infrared measurements that were incompatible with any smooth fit (i.e, the 
measurements indicated minima in the far infrared, which is not a physically plausible behavior), 
suggesting that the issue is the data. Thus, even for the individual systems where we find relatively 
weak star formation, it is an essential part of the SED fits. This conclusion is consistent with our 
finding in stacking the sources not detected individually that a star formation spectral component 
is present on average, although weaker than for the individual Herschel -detected objects. 

A source of systematic error in the decompositions is the probable inclusion of some far infrared 
emission due to star formation in the Elvis AGN template. We have determined a bounding case 
(maximal level of star formation) for this effect as described in Appendix C. Around 160 gm, 
this estimate attributes 75% of the template emission to star formation, so it is impossible to 
apply a substantially larger correction. We have repeated the SED decomposition with this star- 
formation-adjusted template and find for typical cases (where the FIR star formation component 
of the decomposition is substantially stronger than the AGN template) that the upward correction 
in the star-forming luminosity is only ~ 10%. Larger corrections apply for the twelve sources listed 
in the preceding paragraph with relatively weak star formation. Based on the star-formation- 
adjusted Elvis template, the individual corrections to the estimated star forming luminosities for 
these systems are 18, 9, 11, 17, 40, 11, 9, 17, 8, 48, 29, and 11% respectively for galaxies 8, 24, 27, 
38, 48, 49, 61, 64, 69, 71, 77, and 90. For the stacked SEDs, the possible increases in the SFRs are 
25, 10, and 46% respectively, for 0.1 < z < 0.7, 0.7 < z < 1.2, and 1.2 < z < 1.9. Applying these 
corrections would increase the necessity for star-forming templates in fitting the Herschel -detected 
objects and would put the stacked results closer to the ones for the Herschel -detected galaxies and 
emphasize that the non-detected galaxies are, on average, similar but modestly fainter in the FIR. 

There is another important conclusion indicated by Figures [4] and [5] and the SED decompo¬ 
sition. For all members of our sample, the 24 /mi flux density is dominated by emission from the 
AGN, whereas the emission in the Herschel bands is dominated by star formation. Thus, the 
sample selection criteria are unaffected by the level of star formation in the host galaxies. The fact 
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that the the Herschel detection rate does not fall significantly with increasing redshift indicates that 
both the AGN and star forming luminosities in the sample increase with redshift at roughly similar 
rates, that is, the host galaxy star formation must be roughly proportional to AGN luminosity (at 
least that at 24 /jrn). 


3-4-3. Warm Excess 


We found a strong excess above the SED decomposition result from 3 to 60 //m (rest-frame) 
for some sources (see Figure [8] as an example); an additional warm component in addition to the 
star formation and AGN templates is needed to obtain a good fit. A similar excess is also found 
in some z ~ 6 quasar SEDs (Leipski et al. 2013, 2014J). A theoretical model of a parsec-scale 
starburst disk (Thompson et al. 2005; Ballantyne 2008) predicts that a warm component heated 
by star formation would emit strongly in this wavelength range. To introduce a minimum of free 
parameters, we added a component with this specific spectrum (Ballantyne 2008, Figure 7) to 
the SED decomposition template library. There are of course alternative possible origins for this 
emission. Figure [8] shows the comparison in one example of the SED fits before and after adding 
the warm component. The emission from the parsec-scale starburst disk reproduces the hot excess 
very well. The total luminosity from the warm component for this source accounts for 56% of the 
total IR luminosity in this example. 


We judge the fits to require the warm component if the observed 12 /im, 24 /xm and 100 /irn 
fluxes are about twice (or more) the fluxes from the SED decomposition only using AGN, stellar 
and star formation templates. The results are summarized in Table [5] There are eight sources 
that require a warm component to achieve a satisfactory fit, with a contribution to the total IR 
luminosity in the range of 30% to 75%. 


We also tested the influence of a warm component on the conclusions from all of the fits where 
it did not appear to be required; the results are also in Table [5j The column for Lsf.ir shows in 
parenthesis the fractional reduction in the luminosity of the star-forming component if the warm 
component is added to the fit; for example, for source 1, the fit is not improved with a warm 
component and there is no change in Tsf,ir, whereas for source 2 the warm component improves 
the fit and reduces Tsf,ir by 10%. In this latter case, the total luminosity captured by the fit 
is also increased with the warm component, indicating that it accounts for measurements that lie 
above the simpler fit. 


For galaxies that are relatively faint in the far infrared, the introduction of a warm component 
can make the optimum fit ambiguous. For example, of the 12 galaxies with relatively weak FIR 
discussed in the preceding section, three (8, 24, and 64) could be fitted with a substantial warm 
component. Nonetheless, the purely star forming FIR fits are also valid, all having probabilities > 
20% of being satisfactory according to the values of x 2 . Given that the warm component seems to 
be prominent at high redshift and/or high luminosity, and that these galaxies have low far infrared 
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luminosities and are at modest redshift, the star-forming template is the preferred fit. 


In summary, the decomposition inputs are surprisingly simple: 1) a fixed AGN SED from Elvis 


et al. (1994) with adjustable foreground screen reddening; 2) a far infrared SED appropriate for 


local LIRGs; 3) a stellar population component — although we allowed a broad of stellar population 
SEDs, the fits always converged on one appropriate for a relatively old population; and 4) the warm 
IR component. 


3.5. Comparisons with Other Samples 

3.5.1. Comparison of Herschel-detected and Non-detected AGNs in the UV and Optical 


As discussed in Sec. 3.2.1, for the Type-1 AGNs there is no obvious difference between the 


Herschel -detected and Herschel- non-detected populations in the 24 pm flux density distribution, 
nor in the redshift distribution. We now expand that comparison to consider any differences in the 
UV and optical. The SDSS quasars are initially selected through a combination of optical colors 


and confirmed by optical spectroscopy (Richards et al. 2002). Richards et al. (2003) found that 


for them the relative color A (</ — i') is a good indicator of quasar redness for redshifts between 
0.3 and 2.2. The relative color is the difference between the measured color of a given quasar and 
the median colors of quasars at the same redshift: a quasar with large A(</ — i') could either be 
intrinsically red or be reddened by dust. 

There is SDSS coverage of 4.51 deg 2 of our survey area. There are 84 SDSS optically selected 
Type-1 AGNs in this area and 185 AGNs in our 24 pm selected sample. 61 AGNs are selected by 
both samples; 23 SDSS AGNs are not included in our sample due to their 24 pm flux densities 
being below 1 mJy. The plot in the upper panel of Figure [9] shows the colors of the 185 SDSS-and- 
24/rm-detected quasars {(]' — i') (corrected for Galactic extinction) compared with that of SDSS 
quasars in general. We determine the median colors of quasars at a given redshift (the solid line) 


using data from the SDSS DR7 Quasar Catalog (Schneider et al. 2010). They represent a quasar 


population that is optically-bright and not or only slightly affected by dust-reddening. The 24^m- 
detected quasars range from this line to being significantly (~ 1 magitude) redder; this behavior 
is independent of far-IR properties. The K-S test shows that the relative color distributions of 
Herschel -detected and -non-detected AGNs are not statistically distinguishable (P-value = 0.948). 
In other words, Herschel -detected Type-1 AGNs are not significantly redder than Herschel- non- 
detected ones in the optical. This indicates that the dust responsible for the Herschel detections is 
not producing any significant dust extinction along the line of sight toward the AGNs. This does 
not contradict our fits that included reddening of the quasar template; it just indicates that the 
dust responsible does not dominate the far-infrared emission. 
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3.5.2. SDSS optically selected Type-1 AGNs 


In Figure 10 we compare the redshift and i' band magnitude distributions for the SDSS and 
24 fim samples. Basically the two samples select sources in the same redshift range. The 24 gm 
sample includes more sources faint in the i! band than the SDSS sample. The 23 SDSS AGNs that 
are not included in the 24 gm sample are evenly distributed over z ~ 0 — 4, and most of them lie 
at the faint end of the i! band magnitude distribution. The majority of these 23 SDSS AGNs are 
detected at 24 /mi, with flux densities in the range of 0.2 — 1.0 mJy. Therefore the bright 24 fin i 
source selection has a large overlap with SDSS selection, but in addition finds many more (by a 
factor ~ 2) Type-1 AGNs based on their optical spectra. 


In the upper panel of Figure 11 we plot the relative color A (g' — i') (see Sec. 3.5.1) as a 


function of redshift for these two samples. Most SDSS AGNs are scattered around the relative 
color A (g' — i') = 0.0. This indicates that the SDSS AGNs in our survey fields are typical of SDSS 
AGNs in general since the the median value of the color g' — i! is calculated from the large SDSS 
Type-1 AGN sample. Our 24 fin i selected sample includes additional red sources not identified by 
optical colors. At 2 < 1, the red colors in our sample may partly arise from the contribution of the 
stellar component in the optical. Above z ~ 1, the 24 fim selection picks up more luminous AGNs, 
and the SED is dominated by the AGN from the UV to the MIR; these red sources probably either 
have strong dust reddening or intrinsically red AGN continua. 


The lower panel of Figure 11 shows the observed-frame [24 fim/i'] flux ratio as a function of 


redshift for these two samples. The 23 SDSS AGNs that our 24 /an selection missed show small 
[24 fim/i'] flux ratios due to their low 24 fin r flux densities, and most of them are above z ~ 1. 
The contribution of the stellar component is also reflected in the [24 fim/i'\ flux ratio for sources at 
z < 1 for both samples. For most of the AGNs at z < 1, the [24 fim/i'] ratios are lower than that 
predicted by a quasar template. At z > 1, the SDSS-selected AGNs are more consistent with the 
quasar-template-predicted ratios with Ay = 0, while about half of the 24 /un-selected AGNs are 
more consistent with the template predicted ratios with Ay = 0.5. The redder color of [24 fim/i'] 
for 24 /^m-selected AGNs may partly be due to dust-reddening or to intrinsically red continua in 
the optical, or partly due to the warm excess in the MIR enhancing the 24 fim flux density. In 
any case, these results demonstrate that 24 fim selection yields a substantial number of red quasars 
that are absent in purely optical selection. 


3.5.3. X-ray Selected AGNs 


We estimated the equivalent AGN X-ray flux to our 1 mJy selection threshold at 24 fan using 
the bolometric conversion from Lusso et al. (2012), obtaining a flux of 10” 14 erg cm -2 s" 1 in the 
0.5-2 keV band. From the number of X-ray sources at a similar detection limit found by Cardamone 
et al. (2008), we expect a total of ~ 260 X-ray selected AGN above this flux level in our total 
field. Our sample includes 205 infrared-selected sources, or corrected for the ~ 18% incompleteness, 
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about 240. However, from Cardamone et al. (2008), we expect the two samples to have different 
properties, despite the near-coincidence in numbers. 


The intrinsic X-ray luminosity in the [2-10] keV band for a typical AGN at z ~ 1 in our 
sample is 2 x 10 44 erg s , converted from a bolometric luminosity (~ 5 x 10 45 erg s' 1 ) using 
the bolometric to X-ray luminosity correction in Figure 9 of Lusso et al. (2012). The failure to 
detect ~ 80% of the sample in the X-ray suggests that some members are moderately absorbed, 
consistent with their selection in the infrared. From these arguments, most of the AGNs in our 
sample are intrinsically more luminous in the X-ray, compared with the X-ray selected, moderate- 
luminosity (Lx = 10 42 — 10 44 erg s _1 ) AGN sample in Mullaney et al. (2012). AGNs in both 
samples have comparable IR star formation luminosities and (specific) SFRs, and all reside in 
massive, main-sequence star-forming galaxies. However, in the optical and NIR, the SEDs of those 
X-ray selected, moderate-luminosity AGNs are dominated by stellar emission (See the average SED 
in Figure 12 of Mullaney et al. (2012)), while for our sample, emission from the AGNs is dominant 
in the optical and NIR, except for some sources in the lowest redshift bins. From Cardamone et 
al. (2008), the X-ray sample is expected to include nearly 50% of sources where the near infrared 
is dominated by stellar emission, whereas the IR-selected sample is expected to be dominated by 
power-law sources in the near IR (See Donley et al. 2008). That is, a significantly higher fraction 
of the AGN luminosity emerges in the infrared for our infrared-selected sample than is the case for 
X-ray-selected samples. 


4. 24 ^m-Selected Type-2 AGNs 

We now describe the Type-2 AGN identified from the same dataset. Many members of the 
Type-2 sample are at relatively low redshift and their AGNs tend to be of lower luminosity than 
the Type-1 objects. After finding the Type-2 objcts, we will identify a subsample that is directly 
comparable with the Type-1 AGNs in terms of redshift, black hole mass, and accretion rate. 


4.1. Type-2 AGN Identification 

The sample of LIRAS Type-2 AGNs is constructed based on the following selection criteria: 

1. Spitzer /MIPS 24 /xrn flux densities above 1 mJy. 

2. Optical spectra showing narrow permitted emission lines (full width at half maximum (FWHM) 
< 1200 km s' 1 ) with high-ionization line ratios. 

3. If z > 0.34, [Nev] AA3347,3427 detected or FWHM([Oiii]) > 400 km s _1 . 

Because the primary identification is based on emission line strengths rather than widths, well- 
calibrated spectra are required. To maintain consistency in the classification, we only used Hec- 
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tospec data; as a result, a few AGN identified in SDSS that were not targeted with Hectospec 
may have been omitted from the Type-2 sample (see Tables 6 & 7). We exclude any AGNs in the 
clusters according to the source redshift. 

A Hectospec fiber diameter subtends l".5 on the sky. At z = 0.3, 1".5 subtends about 6.7 kpc 
and at z = 0.6, 10 kpc, so the Hectospec fiber includes substantial light from the host galaxy. The 
AGN emission lines can therefore be contaminated by stellar absorption lines from the galaxy. Fol¬ 
lowing Hao et al. (2005), we used the following procedures to subtract the host galaxy contribution 
before measuring the AGN emission lines. First, we select a sample of 212 high S/N spectra of 
pure absorption-line galaxies from SDSS Data Release ?[^} Second, we apply principal component 
analysis (PCA) to construct a library of galaxy absorption-line spectral templates. Third, we fit a 
galaxy template, an A-type star template to account for the young stellar population in the host 
galaxy, and a power-law component proportional to A -Q for the nonthermal component from the 
AGN. A x 2 minimizing algorithm was used to determine the synthetic stellar absorption spectrum. 
Only after stellar and power-law continuum subtraction from all the spectra do we measure the 
emission lines. 


We fitted the following emission lines for each spectrum: Ha, [N n] AA6584, 6548, H/3, [O ill] A5007. 
We rejected all objects with broad components (FWHM > 1200 km s _1 ) in their emission lines 
(i.e., in Ha, H/3, or Mgn 2800). The minimum [OIII] line width criterion of 400 km s -1 was based 
on the fitted width with no allowance for the spectral resolution. Given the resolution of R ~ 1000, 
it corresponds to a threshold of about 270 km s -1 for the intrinsic quasar line width. It should 
therefore not eliminate legitimate AGNs (Brothert on]l996) but protects against inclusion of chance 
anomalous star forming galaxies (e.g., Stanway et al. 2014). We also rejected weak emission-line 


galaxies, i.e., the equivalent width of one of Ha, [Om], or H/3 was required to be greater than 3 

A. 


There are several line flux ratio criteria to distinguish Type 2 AGNs from other narrow emission 
line objects (e.g., Kewley et al. 2001; Kauffmann et al. 2003). Here we use the one from Kewley 
et al. (2001) for objects at z < 0.34, 


/[Om]A5007\_061_ 

6 V H/3 ) log([Nil]/Ha) - 0.47 ' 


( 2 ) 


Since [N n] is redshifted out of Hectospec spectroscopic range at z > 0.34, we use the following 
(Zakamska et al. 2003) for objects at 0.34 < z < 0.76: 


, / [O ml A5007 \ 

log {—m—) 


> 0.3. 


( 3 ) 
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A few AGN at z < 0.34 are identified from the BPT diagram even if their [O hi] lines are narrower 
than 400 km s _1 . Selection of Type 2 AGN with 2 > 0.76 is not possible with our spectra since 
[O hi] is redshifted out of the spectroscopic range. 


The upper panel of Figure [12] shows the emission-line diagnostic diagram for Type-2 AGNs 
at z < 0.34 selected in our sample using Equation [2] The lower panel of Figure 12 shows the 
distribution of the [OIII] to H/3 line ratio for all selected Type-2 AGNs at 0.34 < 2 < 0.76. 


4-1-1. Results 


We identified a total of 85 24 //in-selected Type-2 AGNs over the 3.6 deg 2 survey area; 55 are 
securely detected at least in two Herschel bands, as listed in Table 6. Figure [13] shows two typical 
SEDs for Herschel- detected AGNs. The remaining 30 sources not detected with Herschel are listed 
in Table 7. The redshifts and key derived parameters of the Herschel- detected objects can be found 
in Table 8. We stacked the signals for Herschel- non-detected Type-2 AGNs in two redshift bins: 
0.0-0.4, and 0.4-0.8. There are 14 and 13 sources in these two bins respectively, after rejecting those 


contaminated by close bright objects. The stacked SEDs are shown in Figure 14 


4.2. Morphologies 


We visually examined the Suprime-Cam images (shown in Appendix A) to classify the Herschel- 
detected host galaxy morphologies (51 out of 55 sources have images of adequate quality). Although 
the images are not of sufficient resolution for a definitive determination of the morphologies of the 
entire Type-2 sample, they allow us to make plausible assignments for most members (summarized 
in Table 8). Twelve of the 55 Herschel- detected galaxies either have no useful imaging data (4) or are 
sufficiently compact that no further morphological information can be derived. Most of the rest (38) 
are early-type spirals, lenticular galaxies, or elliptical galaxies. Only five are probable interacting 
systems, although a few of the ellipticals and lenticular galaxies also show hints of distortion and 
interaction. Therefore the majority of the AGNs reside in normal-appearing spheroidal and bulge- 
dominated galaxies. This result is consistent with the results of Povic et al. (2012) for a sample 


of X-ray selected AGN at z < 2.0 and with those reported by Villforth et al. 
CANDELS fields at z ~ 0.7. 


(2014) for the 


4.3. Type 2 SED Decomposition 

The Type-2 AGN sample does not extend to z > 0.8 because the critical emission lines move 
outside the range of our spectra; in fact, from the trend of detections with redshift, the sample 
becomes progressively less complete above z = 0.6. The sample includes many more galaxies at low 
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redshift (z < 0.3) than for the Type-1 sample. Many of these low redshift galaxies appear to be 
dominated by star formation at 24 ^m, with AGNs both of relatively low luminosity (because of 
the low redshift) and, by themselves, not as bright as 1 mJy at 24 gm. To make these statements 
quantitative, we need to carry out the decomposition of the spectral energy distributions. Only 
then can we determine which sources can be compared directly with the members of the Type-1 
sample. 


4-3.1. Decomposition Approach 


For both the individual Herschel -detected sources and the stacked results, we use SED de¬ 
composition to disentangle the AGN and star formation contribution in the FIR. Based on the 
arguments in Sections 3.2.3 and 3.3, we assume that star formation dominates the signals in the 
Herschel bands. Specifically, we model the observed SED as a linear sum of a stellar component, 
a star formation component, and an AGN component. We use Levenberg-Marquardt least-squares 
fitting to find the best stellar, star formation, and AGN templates and their normalizations. 


The stellar population and FIR star-forming galaxy templates were identical to those used 
with the Type-1 sources. The GALEX data show UV excess emission in the majority of the Type 
2 AGNs that cannot be produced by an old stellar population. This enhanced blue color is also 
reported in Sanchez et al. (2004) for early-type AGN hosts at 0.5 < z < 1.1. SDSS images of 
Seyfert 2 galaxies at z < 0.2 show UV emission from young stars in the outer regions of the host 
galaxies (Kauffmann et al. 12007). Because the stellar population SEDs did not allow for two distinct 
episodes of star formation widely separated in time, we added SEDs for a second population of very 
young UV-bright galaxies of age 0.1 Gyr to 1.0 Gyr to fit the UV emission. 


The strong extinction associated with Type-2 nuclei makes it impossible to use a single template 
to fit their SEDs. We employ the numerical AGN templates from Fritz et al. (2006), which include 
cases with heavy absorption. These models assume a central point-like energy source with a broken 
power-law SED surrounded by a smooth dust distribution, and then solve the radiative transfer 
equation. Templates generated using this method depend on the dust distribution geometry and 
composition, and the inclination of the torus toward the observer. We put constraints on the 
AGN template library based on observations of the Seyfert 2 galaxy silicate 9.7 /jrn absorption 
features, for which it is found that ( Ff — F c ) / F c > —0.85, where Ff and F c are the observed flux 
density and underlying continuum flux density at the minimum of the silicate absorption feature, 
respectively (Shi et al. 2006). Therefore, we do not use AGN templates with silicate 9.7 fim 
absorption (Ff — F c )/F c < —0.85. For comparison, Nenkova et al. (2008) have calculated models 


for clumpy dust distributions; the comparison of smooth and clumpy dust models in Feltre et al. 


(2012) shows that, although the two types of model give different outputs, for our purposes they 
are equivalent. For example, the Nenkova et al. (2008) models limit the silicate absorption depth, 
which we did also by imposing an additional constraint on the Fritz et al. (2006) ones. 
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Samples of the deconvolutions are illustrated in Figures 13 and 14. Table 8 lists the derived 
parameters for all 55 Herschel- detected sources. 


4-3.2. Relative Roles of Star Formation and AGN at 24 /ira 


All of the Type-1 AGNs are dominated by emission by the AGN at 24 /am (including the 
four at z < 0.3^ We classify a Type-2 AGN as AGN-dominated if the flux arising from AGN 
component at 24 /jin is larger than that from the SF component; otherwise, it is defined as 24 fin r 
SF-dominated. Figure [13] shows examples of these two classifications. There are 17 AGN-dominated 
and Herschel- detected Type-2 galaxies, about 1/3 of the Herschel- detected sample. There is an 
approximate divide in this behavior at z ~ 0.3. Above this redshift, there are 50 Type-2 galaxies in 
our sample, of which 37 (74%) are dominated (or tied) by emission by the AGN over that from star 
formation at 24 fin i (We include the Herschel- nondetected galaxies in this sample, since normalizing 
a star forming template to the Herschel upper limits shows all of these to be AGN-dominated). 
However, for the 35 cases at z < 0.3, only 13 (37%) are AGN-dominated at 24 fim (including 
the Herschel- nondetected cases). In Section 3 we showed that selection at 24 pan yielded a large 
number of Type-1 AGNs; it appears that for z > 0.3, selection at this wavelength is useful to 


generate candidate lists that are relatively unbiased in terms of AGN type (see also Mateos et al 


(2013)). 


Because the 24 fan selection works relatively well at z > 0.3 in finding the most luminous 
AGN, we use it to compare the incidence of Type-1 and Type-2 sources. There are 50 Type-1 
AGN with 0.3 < z < 0.8 in the 5.2 square degrees surveyed for them; normalizing by surveyed 
area, we expect 35 in the 3.6 square degrees surveyed for Type-2 sources. In fact, we have found 
37 dominant Type-2 sources. That is, the numbers of Type-1 and Type-2 quasars in this redshift 
range are similar. This result confirms the conclusion of Reyes et al. (2008]), but with the initial 
selection on a completely different basis than the extinction-corrected [OIII] luminosity used in that 
work. 


4.4. Definition of the High Luminosity and Comparison Samples 

4.4H. Sample Definition 

We now derive a subsample of Type-2 objects suitable for comparison with the Type-1 AGN. 
The Type-1 sources are very luminous, with massive black holes (77/91 = 85% have M bh > 1 X 10 8 


9 In some cases, the warm component is dominant. The origin of this component is not known, but it appears to 
be unique to AGN. 
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Mq) and accreting at rates close to Eddington (66/70 = 94% at > 3% of the Eddington rate] 10 
To compare with them, we need to define a suitable sample of the Seyfert 2 AGN, namely those 
indicated to have Mm > lx 10 8 M 0 and that are accreting at a minimum of 3% of the Eddington 
rate, leading to a minimum bolometric AGN luminosity of 10 11 Lq. We describe these objects as 
the high luminosity sample (HLS) as indicated in Table 8. The HLS consists of 17 objects, all but 
two at z > 0.3. All but four of the 107 Herschel -detected type 1 galaxies are also at z > 0.3. Yan 


et al. (2013) show that the incidence of star-forming galaxies bright enough to be within our 24 


/jin selection is low at z > 0.3, simplifying the task of identifying AGN. Therefore, for the primary 
Comparison Sample with the Type-1 objects, we require z > 0.3; the 15 members of this sample 
are also flagged in Table 8. 


Not surprisingly since both metrics emphasize high AGN luminosity, 13/17 of the AGN- 
dominated sources also belong to the HLS. By definition, the 15 sources in the Comparison Sample 
are all members of the HLS, but the Comparison Sample also includes 14/17 of the AGN-dominated 
examples. Thus, the various methods for isolating the most luminous AGNs largely overlap. How¬ 
ever, because its membership is not linked to the host SFR and its threshold AGN IR luminosity 
is matched to that in the Type-1 sample, the Comparison Sample is best suited to complement the 
Type-1 sample. 


Possible Biases in the Comparison Sample 


The members of the Comparison Sample virtually all fall in the range where we identified 
the AGN by the ratio of [O hi] and H/3 line fluxes. We now consider the reliability of this iden¬ 
tification procedure. Figure 1151 shows the correlation between the line luminosity, ^o m j> and 
the AGN total luminosity, the latter from our SED decomposition. The two luminosities corre¬ 
late as Lagn oc -^Q 4 m j j even though it is generally believed that the strength of the [O ill] line 
should be proportional to AGN luminosity. At the higher end of the redshift range, sources have 
higher L^Q m j than pure proportionality predicts. Since LjQ in j traces ionizing photons that can 
be created by star formation as well as AGNs, one possible reason for the L|-Q m j excess at higher 
redshift is that the field of view of the fiber of the spectrograph includes significant amounts of 
[O hi] emission from star formation in the host galaxy, as discussed further in Xu et al. (2015). 


The possibility, particularly at high redshifts, that our [O hi] measurements are contaminated 
by the host galaxies could result in a bias against AGNs in host galaxies with very strong star 
formation, since they might be expected to have reduced ratios of [O hi] to H/3 and thus miss our 
selection criteria. However, we believe this is not a problem for a number of reasons. First, we have 
searched for candidate contaminated systems at z > 0.34 through the entire spectroscopic sample, 


10 The denominators for both of these percentages are based on the number of objects with suitable measurements 
and do not represent the entire sample. 
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by identifying galaxies with [Om] FWHM > 400 km s^ 1 , 1.5 < [Oiii]/H/3 < 2 (corresponding 
to 0.176 < log([0 m]/H/3) < 0.3, below the selection threshold in Equation [3]) , and 24 /mi flux 
density > 1 mJy. Because all of our candidate galaxies have masses > 3 x 10 10 Mq, this selection 
procedure would have identified all Type-2 AGN in our stellar mass range that would satisfy the 
MeX criteria (see Figure 4c in (Juneau et al. 2011)). We found only one candidate. This low 


yield is consistent with our only identifying 5 out of our sample of Type-2 AGN above z = 0.3 
with [O hi]/H/ 5 between 2 and 3; it appears that our initial 24 /mi selection generally yields AGN 
with relatively large [O m]/H/3. This result suggests that there are very few candidates that might 
have missed identification as Type-2 AGN because of contamination. The low yield with a relaxed 
log([0 iii]/H/3) threshold also indicates that the Type-2 sample is nearly complete to an AGN flux 
density of 1 mJy, at least for z > 0.3. 

Second, consistent with this conclusion, if contamination were introducing significant biases, 
one would expect that the Herschel -detected systems would have a tendency to have low values of 
[O III] /H/3 because they have relatively strong star formation, but the lower panel of Figure 12 does 
not show a strong effect. 

Third, it appears that our AGN samples include all potential contaminated galaxies. We have 
determined that the ratio of [O ill] to 24 /mi flux density is roughly the same or slightly higher 
for AGN compared with star forming galaxie^j In addition, the value of [Oiii]/H/3 intrinsic to 
AGN is often significantly higher than our adopted theshold of 2. As an example, for the sample 

(2010]), the median ratio is 9, while the sample of Juneau et al. 


compiled by LaMassa et al. 


(2011) has a typical ratio of 4. Taken together these results show that any host galaxy containing 


an AGN with an intrinsic flux density > lmJy at 24 /mi (i.e., above the luminosity threshold for 
our AGN samples) plus star formation sufficiently vigorous to contaminate the [O III] /H/3 ratio 
enough to cause it to fall below our threshold would have a total signal at 24 /mi well above 1 mJy. 
However, above 2 mJy, the AGN in our samples count for all of the detections, leaving no room for 
a population of luminous AGN in very luminous star-forming galaxies. 

Figure[l5|is not the first finding of a departure from the expected 1:1 relation between [O hi] and 
bolometric AGN luminosity in the direction of an increasing [OIII] lluminosity for more luminous 


AGN. LaMassa et al. (2010) found similar behavior relative to 12 /mi luminosities; Shao et al. 


1 For AGN, we used the sample from 


Diamond-Stanic et al. 


(2009) (nuclear 24 /an flux densities were provided 


by A. Diamond-Stanic, private communication). We found an average value for the flux in [OIII] (in aW) over the 
flux density at 24 fim (in Jy) of 513 ± 80 (492 ± 69 for type 2 AGNs and 548 ± 118 for type 1). For star forming 
galaxies, we utilized the MIPS 24 /im measurements and the ’’radial strip” line results for the SINGS sample ([Dale 


et al. J2005 


Moustakas et al. ~|[2010) to find a ratio of 457 ± 114 and the integrated galaxy spectra from Moustakas 


& Kennicutt | (2006) with IRAS 25 /im data to obtain 426 ± 15. We have also used the ’’radial strip” spectra from 


Moustakas et al. 


(2010) for galaxies with My < -20 to find an average value of [O ill]/H/3 = 0.98 for luminous 
star-forming galaxies. Relaxing the luminosity threshold to My < -19 has little effect: the average is then 1.03, so 
the value is not strongly sensitive to galaxy luminosity (and the accompanying range of relevant metallicity, based on 


the luminosity-metallicity relation). Caputi et al. (2008) find a similar average ratio, while the work of Moustakas 


& Kennicutt (2006) yields a value of 0.82, again in good agreement. 
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(2013) also saw this behavior when comparing the [OIII] and 22 /j,m luminosities for a large sample 
of AGN from the Sloan Digital Sky Survey; and Hainline et al. (2013) report a similar departure 
from a 1:1 relation using 8 ^m luminosity as an indicator of AGN luminosity. However, Shao et 
al. (2013) show a 1:1 relation between [OIII] and 4.6 /jm luminosity. Taken together these results 
imply discrepancies in measuring the AGN luminosities from different single-color infrared bands. 
It is therefore of interest that we find the effect based on the bolometric AGN luminosity, rather 
than an estimate of the luminosity based on a single spectral band. 


5. Intrinsic Properties 

With the SED decompositions, along with other derived properties of the sources (e.g., line 
widths), we now derive the basic physical parameters of the sources in our sample. 


5.1. IR Luminosities and Star Formation Rates 


The SED decomposition disentangles the contributions from different source components. The 
infrared luminosity from the star formation component (Asf,ir) is integrated over the rest-frame 
8 — 1000 /im range of the best-fit star formation template. The infrared luminosity from the 
AGN component (Aagn,ir) is integrated over the same range of the rescaled AGN tempi at (j^j By 
integrating the full Elvis template, we set the total AGN luminosity of the Type-1 objects to be 
5.28 times their infrared luminosities. The Type-2 AGN bolometric luminosities are taken to be 
the rescaled total intrinsic luminosity of the best-fit AGN template for each source 13 The total 
infrared luminosity (Air) is the sum of Lsf.ir. and Aagn.ir- The star formation fraction (Fsf) 
is defined as Asf,ir/Air. Of the Herschel -detected Type-1 sources, 21% have Fsf > 75%, 47% 
have 50% < Fsf < 75%, and 32% have F$f < 50%. The corresponding values for the Type-2 
Comparison Sample are 60 ± 20%, 27 ± 13%, and 13 ± 10%, similar within the poor statistical 
weights of the latter (particularly allowing for the lower typical redshifts of the Type-2 objects). 
The star formation fractions of the stacked SEDs of Herschel- non-detected Type-1 AGN (in three 
discrete redshift bins: z =0.1-0.7, 0.7- 1.2, and 1.2-1.9) are all below 40%. 


“The calculated luminosities are uncertain for a number of reasons, such as: 1.) the underlying assumption that 


the emission by the central engine is isotropic, despite its complex geometry and optical depth (e.g., Koratkar & Baes 


(19991; 2.) the contamination of the AGN template in the FIR by star formation (see Appendix C); and 3.) the 
inclusion of both the optical-UV-X-ray and the infrared components of the SED ( Marconi et al.|2004 |. It is difficult 
to make quantitative estimates of these effects, but other than the first, they appear to be modest (i.e., < a factor of 
two for the second two together (see 


Marconi et al. 


(2004) and Appendix C). 


^Uncertainties for them include: 1.) the differences between clumpy and smooth models (Feltre et al. 2012); 2.) 
in our fitting, the torus opening angle is poorly constrained; 3.) variability; and 4.) the underlying assumption that 
the emission by the central engine is isotropic, despite its complex geometry and optical depth (e.g., Koratkar & 
Blaes (1999)) 
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The star formation rates are calculated from LgF,iR, using the relation in Kennicutt (1998), 
adjusted for a “diet” Salpeter initial mass function (IMF) (Bell et al. 2003) from the original 
Salpeter IMF, i.e., 


SFR 
M 0 yr- 1 


= 1.2 x 1(T 10 


/ TsF,IR 

V L0 


( 4 ) 


The adopted IMF reduces the proportion of low mass stars to resemble, for example, the Kroupa 
IMF, and puts the SFRs on the same scale as our mass estimates in the following section. Tsf,ir> 
ranges from ~ 10 10 to 3 x 10 12 L 0 for the Herschel -detected galaxies; the average value for the 
stacked SEDs of the non-detected galaxies are several times lower. Nonetheless, star formation 
activity must be common even for the AGN hosts not individually detected by Herschel. However, 
as for the local sample of PG quasars (Shi et al. 2014), it is possible that there are a number 
of quiescent galaxies among those we stacked, and therefore that elevated star formation is not 
ubiquitous. 


5.2. Virial Black Hole Masses and Eddington Ratios 


Type-1 AGN black hole masses, M., have been measured directly by reverberation mapping 
(Blandford & McKee 1982; Peterson 1993; Kaspi 2000), but it takes years to obtain results using this 
technique. However, reverberation mapping has also provided empirical scaling relations allowing 
us to estimate black hole virial masses efficiently from the quasar continuum luminosity and broad 
emission line widths, e.g., H (3 (4861 A), Mg II (2800 A) , and C IV (1549 A). We used the moderate 
resolution (~ 6 A, corresponding to 300-400 km s _1 ) Hectospec spectra to determine the FWHM 


of the broad emission lines. We followed the procedures in Vestergaard & Wilkes (2001) (for Mg 
il), and Peterson et al. (2004]) (for H (3 and and Civ) to fit these lines and measure the FWHM of 


the broad component. We took the mass-scaling relationship from Vestergaard & Peterson (2006) 


(for H/3 and Civ) and from Vestergaard & Osmer (2009) (for Mg il) to estimate black hole masses. 


In Appendix [Cj we list the three mass-scaling relationships we used, and show three examples of 
fitting results for Mgn, Civ, and H/3, respectively^] The measured FWHMs and estimated BH 
masses are listed in Table[5j 


The Eddington Luminosity from a black hole with mass M. powered by spherical accretion is 

( 5 ) 


Le = ( 4TGcm - | M.. 


V Ce 

We obtained the AGN total luminosity from the SED decomposition, and calculated the ratio of 
AGN luminosity to Eddington luminosity. Of the Type-1 AGN, 94% emit at <13% of the Edding- 


14 If both H/3 and Mgll are available, we adopt H/3; and if both Mgn and Civ are available, we adopt Mgn. 
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ton Limit. The distribution of the bolometric luminosity as a fraction of the Eddington limit is 
consistent with that of the SDSS quasars (jMcLure & Dunlop] 2004]). 

Assuming the local stellar mass (M*) and black hole mass (M.) correlation (i.e., M * ~ 700M.; 
e.g., Bennert et al. 2011; Cisternas et al. 2011; Scott etal. 2013), we calculate the black hole 
masses, Eddington luminosities, and Eddington ratios of the Type-2 AGNs, with results shown 
in Figure 16 The 24 fim SF-dominated Type-2 AGNs have slightly lower ratios than the AGN- 


dominated ones at all redshifts. The 24 //rri AGN-dominated Type-2 galaxies emit close to 10% of 
the Eddington rate (14/17, or 82% emit at <13% of the Eddington rate); for z > 0.3, their behavior 
is similar to that of the Type-1 AGNs. Therefore, as expected (e.g., by the unified model), the 
behavior of the nuclei of the 24 /im AGN-dominated Type-2 galaxies is consistent with that of the 
Type-1 sample. The lower Eddington ratios in the star-formation-dominated galaxies are expected, 
given that they have not been selected strictly on AGN luminosity. 


5.3. Stellar Masses of AGN Host Galaxies 


5.3.1. Stellar masses from SEDs 

Based on the SED decomposition, we can estimate the host galaxy stellar masses. Because the 
details of the stellar spectrum are difficult to disentangle from the AGN emission, we base the mass 


estimate on the NIR stellar luminosity, which has been shown to be an accurate approach (e.g., Me- 
Gough &; Schombert 2014). We use the relation between stellar mass, M*, and K-band luminosity, 


Lfc, for local field galaxies (Bell et al. 2003): (logio(M*/Lk) = —0.42+0.033 logio(M c h 2 /MQ), where 
M c h 2 is 10.63 averaged over all galaxy types (10.61 for early-type galaxies and 10.48 for late-type 
galaxies). The masses assume a ’’diet” IMF, defined by Bell et al. (2003). We need to be sure that 


our mass estimates are on a consistent scale with other approaches. This would be straightforward 
if the host galaxies were normal early-types, but many of them have anomalously blue colors (e.g., 


Floyd et al. 2013). We therefore compare with a wide range of masses that include galaxies with 


a range of colors. We find that the masses are consistent with those using SDSS KCORRECT 
(Blanton Sz Roweis 2007 15 for local galaxies (private communication, Krystal Tyler; See Figure 


Although photometrically determined masses can be subject to significant systematic errors, 
the agreement on average between our approach and the masses derived from full photometry puts 


our masses on the same scale as, e.g., those of Elbaz et al. (2011) and allows direct comparison of 


the host galaxy behavior in our sample with the held galaxy behavior described in that paper. The 
stellar component in the Type-2 galaxies is more accessible to our fitting than for the Type-1 cases, 
and we could attempt more detailed models. However, for consistency in comparing the samples, 
we use the same approach. We obtain the estimates of stellar mass in Table 8. Our AGNs reside 


15 Also see http://howdy.physics.nyu.edu/index.php/Kcorrect The SDSS KCORRECT stellar mass is based on 
the Bruzual-Charlot stellar evolution synthesis and makes use of the multi-band SDSS photometry. 
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in very massive galaxies with stellar masses around 10 11 Mg. 

As a stellar population ages, its luminosity declines as its more massive stars die; i.e., a fixed 
K-band luminosity corresponds to smaller stellar mass at higher redshift. This passive evolution 
must be accounted for in estimating the masses of the stellar populations in high redshift galaxies 
(e.g., Drory et al. 2003, 2004; van der Wei et al. 2006; also see van Dokkum & Franx, 2001). We 
assume the AGN host galaxies evolve passively and follow van der Wei et al. (2006) to correct for 
this systematic evolution of the host stellar luminosity (i.e., A In (M*/Lk) = (—1.18 ± 0.10)z). 
This correction can be applied to galaxies from the local epoch to z ~ 1.2, where the correction 
factor is equal to 4.1. For z ~ 1.2 — 2, we keep the value of the correction factor at 4.1. 


Since the rest-frame J-band is near both the peak of the stellar SED and a minimum of the 
AGN SED, we use this band to quantify the stellar component output. Our fits constrain the J-band 
flux from the stars well for most low-z sources, where the AGNs are of relatively low luminosity 
and do not dominate in the rest NIR. At higher redshifts, we can usually only obtain upper limits 
for the stellar fluxes. We ran a simulation to test to what level we can trust the stellar flux from 
the SED decomposition. First, we renormalized the stellar and AGN SED templates to the desired 
flux ratio in the rest-frame J-band. Second, we applied dust extinction selected randomly over the 
range Ay = 0 — 1.0 to the AGN template and then added the two templates. Third, we convolved 
the bandpass transmission curves with the combined templates to simulate the photometry that 
we used for the SED decomposition. We added random noise to the simulated photometry in 
all bands, assuming a standard deviation of 20% in consideration of the photometry errors, AGN 
variability, and that the data in different bands were probably taken in different years. Fourth, we 
ran the SED decomposition procedures on this simulated photometry and calculated the recovered 
flux ratio of the AGN and stellar components in the rest-frame J-band. The input flux ratio was 
set to six discrete values: fluxAGN,j/fluxst e ii ar ,j = 0-5,1,2,3,4, 5, and the calculation was repeated 
10,000 times for each value. We compare the input and recovered flux ratios in Figure 18 If the 


stellar flux is twice the AGN flux in J band, 99% of the sources can be recovered accurately. If 
the stellar flux is equal to the AGN flux in J band, there is a larger scatter of the recovered flux 
ratio, and on average, the stellar mass is overestimated by about 20%. However, > 95% of the 
sources are recovered within a factor of two. If the stellar flux is below the AGN flux in J band, the 
errors in the stellar flux are large. Based on this result, if the rest-frame J-band flux of the stellar 
component is equal to or above that of the AGN component, we can compute a valid stellar mass. 
If the J-band flux of the stellar component is smaller than that of the AGN component, we use the 
J-band flux of the AGN component to assign an upper limit to the mass of the stellar component. 


The J — K colors of early-type galaxy stellar populations are very similar, so the rest-frame 
K-band luminosity can be taken to be 0.85 times the J-band luminosity. Therefore we use 0.85 
times the stellar component J-band flux to compute stellar masses, or of the AGN component to 
estimate stellar mass upper limits. We then obtain the estimates of AGN host stellar mass as 
tabulated in Table 5. 
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5.3.2. Indirect determination of stellar masses 

We now estimate stellar masses from the black hole-stellar bulge relation. These estimates 
allow us to 


1. extend the study of AGN host galaxies to a significant number at z > 1 

2. test the passive evolution assumed to correct our mass estimates from observed near infrared 
fluxes 

3. investigate the possible bias toward massive host galaxies because requiring them to be suf¬ 
ficiently bright in the NIR to outshine the AGNs for photometric mass estimation will favor 
ones with massive hosts, at least at high redshifts 


To lay the foundation for indirect mass estimates, we: 1.) examine the possible extent of evolution 
of M./M* over the relevant redshift range, 0 < z < 1.8; and 2.) calibrate the masses derived from 
M. against those from near infrared luminosity obtained in the preceding section. These two steps 
let us determine the maximum plausible deviations of the derived stellar masses from a nominal 
“best estimate.” 


The great majority of luminous AGN are in galaxies with early-type morphologies (e.g. McLeod 


fe Rieke~||1995 Floyd et al. ||2004 ). For such galaxies, the local value of M./M* is well determined 
for galaxies with M. > 3 X 10 7 M 0 (Kormendy & Ho 2013). The majority of our AGN samples 


with z > 0.3 have M. above this threshold, within the range where M./M* is well behaved. 

Most investigators agree that, within the errors, there is little evolution in the M./M* ratio 


from z = 1 to z = 0 ( 

Peng et al. 

2006a 

Shen et al. 

2008 

Somerville 

2009. 

Cisternas et al. 2011; 

Zhang et al.||2012 

Salviander et al. 2013| ?). A small number of studies suggest some evolution in 


this range but are inconclusive regarding its significance ( Woo et al. ||2008 Canalizo et al. 2012). 
We will assume that the local ratio holds up to z ~1.2. For z between 1 and 2, the indications 


al. 

2010, 

Schulze & Wisotzki 20R 

) to evolution by a factor up to about four (at z = 2) (Peng et 

al. 

2006a; Trakhtenbrot Sz Netzer 

2010, 

Merloni et al. 

2010 

Decarli et al. 2010| Bennert et al. 


2011). Particularly at z > 1, there are selection effects that bias the apparent evolution upward, 


correction for which reduces it significantly, to a factor of two or less (Lauer et al. 2007, Shen & 


Kelly 

2010 

Schulze <fc Wisotzki 

2011, 

Portinari et al. 

2012 

16 


However, these same selection 


effects, e.g. the bias due to luminosity selection toward relatively massive black holes (Lauer et al 


2007), may apply to our use of the M./M* ratio to estimate stellar masses 77 ! so for our sample, 


we will consider the possible extreme value of M./M* at z = 1.8 to be 4 times the local value. 


10 The factor of two is also consistent with the conclusion that half of the stellar mass observed today has formed 
since z = 1.3 (Madau & Dickinson J2014I, assuming that any black hole growth over this period is negligible. 


17 


See 


Matsuoka et al. 


(20141 for an example of about a two times bias for galaxies of similar mass to ours. 
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To calibrate stellar masses from M./M* against masses from near infrared photometry, we 
assume the M./M* relation has no evolution from the local value M* ~ 700M. up to z ~ 1.2 (e.g., 
Bennert et al. 2011; Cisternas et al. 2011; Scott et ah 2013). At 0 < z < 1.2, there are 28 AGNs in 


our sample that have both M. (derived from broad line width; See Section 5.2) and M* (estimated 
via K band luminosity). In Figure 1£^, we compare the K-band derived M*’s and those from 


M./M*. The large scatter is not surprising given the scatter in the cr-M relation (Kormendy & Ho 
2013) and in the mass determinations from photometry (e.g., ( Shapley et ah |[2001 Savaglio et ah 
1|2005 Kannappan &; Gawiser||2007 )), plus the additional uncertainties indicated by our simulation 
of the deconvolution uncertainties in near infrared fluxes. The K-band derived M* is on average 
two times higher than the bulge stellar masses predicted by the local M./M* ratio. This offset is 
redshift independent, which indicates our correction for passive K-band evolution is roughly correct. 

The offset could arise if the galaxies have substantial disks, or if our photometric or BH masses 
have small systematic errors. However, there is also a selection bias toward relatively massive 
galaxies that have sufficiently bright NIR stellar fluxes to outshine their AGN and allow mass 
estimates from their SEDs. Approximating this bias by assuming a normal intrinsic distribution 
with all of the cases with mass estimates from SED fits coming from the upper side indicates an 
offset by a factor of 1.7, in satisfactory agreement with what we find. Thus, the indirect stellar mass 
estimates serve the important function of removing this source of bias from our sample. The scatter 
is 0.38 dex rmo We can estimate the intrinsic scatter in our near-infrared-based mass estimates 


of about 0.15 dex from Figure 17 corrected for the scatter in the masses with full photometric 
fits (?). If anything, 0.15 dex may be a low estimate (Courteau et al. 2014). We need to add 


the scatter due to having to measure the near infrared fluxes from deconvolution of SEDs with 
significant contributions from the AGNs. The resulting total scatter is expected to be at least 0.2 
dex. Subtracting this value quadratically from 0.38 dex, we estimate that the masses determined 
from M./M* will have an intrinsic rms scatter of 0.32 dex. Given the uncertainties we estimate the 
possible errors as 0.3 dex toward low values (from the rms scatter) but 0.6 dex toward high values 
above z = 1.1 (from possible evolution and/or selection effects), relative to the nominal values 
assuming the local M./M* ratio. All the AGNs with stellar masses from photometry reside in very 
massive galaxies with stellar mass around 1 — 4 x 10 11 M 0 . The masses estimated indirectly are 
also generally within this range. 


^One extreme outlier has been omitted; such outliers are also seen in other samples (Kormendy & Ho 


9 This value is plausible given the expected errors in single-epoch black hole mass estimates found by 


2013) 


Vestergaard 


& Peterson (2006), corrected for scatter in the reverberation mapping masses (Onken et al. 20041, or the single 


epoch error limits estimated by Denney et al. (20091 
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6. Summary 

We studied the properties of a sample of 24 /jin-selected, spectroscopically-identified AGN and 
their host galaxies, using a multi-wavelength dataset from Chandra, GALEX, SDSS, UKIRT, 
WISE, Spitzer/ MIPS, and Herschel. Typical luminosities for these AGN are above 10 45 ergs s^ 1 
(~ 2 x 10 11 Lq), and they generally lie between z of 0.3 and 2.5. We use SED decomposition from 
the optical to the FIR to estimate the AGN luminosities, SFRs and stellar masses of the AGN 
hosts. 

We summarize the results from this study as follows. 

1. About 50% (107 out of 205) of the Type-1 AGNs in our sample are individually detected by 
Herschel. Among these AGNs, 68% show high levels of star formation (the star formation 
activity contributes over 50% in the FIR). Herschel non-detected AGNs were studied using 
stacking analysis. On average, they have a similar level of AGN luminosity and similar optical 
colors, but the average star formation activity is several times lower compared with AGNs 
individually detected by Herschel. 

2. Similarly, about 65% (55 out of 85) of the Type-2 AGNs are individually detected by Herschel. 
However, these objects tend to be at relatively low redshift and some of the detections are 
a result of vigorous star formation, not nuclear activity. We have defined a sample of 15 
Type-2 AGN with properties (M bh-, Eddington ratio, and redshift) that make them directly 
comparable with the Type-1 sample. 

3. The FIR-detected Type-1 AGNs and matching Type-2 ones reside in massive galaxies (~ 
1 — 2 x 10 11 Mq). They harbor supermassive black holes of ~ 3 x 10 8 Mq, which accrete at 
~ 10% of the Eddington luminosity. 

4. A warm excess in the MIR was found for eight Type-1 AGNs compared with a local quasar 
template. This warm excess can be prominent at higher redshifts but is not seen in low 
redshift quasars. It is not clear whether it changes due to evolution, or whether the warm 
excess is confined to very luminous quasars. 

5. The 24 yum-selected sample of Type-1 AGNs includes about twice as many objects as are 
identified through the SDSS, including the majority of the SDSS identifications. The ad¬ 
ditional objects have redder optical colors than typical SDSS quasars, due to reddening or 
intrinsically red quasar continua. 

6. As also found, e.g., by Hainline et al. (2013), the strength of the [OIIIJA5007 line increases 
more rapidly than proportionately to bolometric AGN luminosity. At relatively high redshift 
(and hence high AGN luminosity), detection of [OIII] emission from parts of the host galaxy 
within the spectrograph fiber may contribute to this effect. 

These results are discussed further in Xu et al. (2015). 
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Table 1. Cluster fields 
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a Type 2 AGNs are selected using MMT/Hectospec 
spectra in these 21 cluster fields. Because the selection 
depends on emission line ratios, we did not select Type 
2 AGNs in 9 cluster fields where the spectra were not 
flux calibrated due to the lack of standard stars. 
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c Some AGNs are broad absorption line quasars, and their broad emission line FWHM is above 1200 km s , but is not well constrained. Some AGNs have 
low S/N optical spectra, and the emission line fitting is not good enough for black hole mass estimate. 



Table 6. Fluxes of 24 /xm-selected Herschel -detected Type-2 AGN sample 1 
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Table 7. Fluxes of 24 /xm-selected Herschel- non-detected Type-2 AGN sample 
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CN CO ^ LO CO b 00 05 


CNCO^iOCDbOOC50HdCO^iCCONOOC50 
t-H t-h t-h t-h i-h t-H t-h t-h CN CN CN CN CN CN CN CN CN CN CO 


J083229.03+654710.6 128.1209420 65.7862700 0.730 0.009 0.023 0.03 0.05 0.32 1.71 

J024713.00-033745.1 41.8041692 -3.6291907 0.173 0.229 0.302 0.18 0.17 0.84 1.04 

J172057.08+263051.4 260.2378232 26.5142879 0.187 0.862 1.009 0.48 0.44 1.46 3.27 

J145733.19+221450.1 224.3883064 22.2472405 0.532 0.137 0.278 0.20 0.15 0.52 1.97 

J133455.26+410724.7 203.7302636 41.1235413 0.633 0.068 0.131 0.17 0.14 0.31 1.36 
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Table 8. Redshifts and derived parameters for the Herschel -detected Type-2 sources 


# 

Source 

LIRAS 

z 

-^tot,IR 

1O 11 L 0 

Lsb.ir 

1O 11 L 0 

Lagn.ir 

1O 11 L 0 

^AGN, total 

1O 11 L 0 

M* 

IO^Mq 

Type d 

1 

J 101756.90+390528.0 

0.054 

0.10 

0.09 

0.01 

0.03 

0.33 

S 

2 

J 101805.64+385009.7 

0.067 

0.21 

0.21 

0.00 

0.03 

0.77 

S0/E 

3 

J 145753.24+222422.7 

0.109 

0.13 

0.13 

0.00 

0.01 

0.76 

C 

4 

J133250.99+501816.1 

0.110 

0.12 

0.11 

0.01 

0.32 

0.31 

S0/E 

5 

J083713.49+150037.4 

0.141 

0.23 

0.22 

0.01 

0.13 

0.80 

E 

6 

J080102.09+355132.2 b 

0.160 

0.34 

0.19 

0.15 

2.14 

0.60 

E 

7 

J 133644.33+405854.9 

0.169 

0.30 

0.23 

0.08 

0.32 

1.00 

S 

8 

J101623.97+385840.1 

0.169 

4.75 

4.72 

0.03 

0.70 

0.68 

E 

9 

J172109.90+263455.1 

0.170 

0.67 

0.60 

0.06 

0.78 

0.55 

I 

10 

J172022.13+263626.6 b 

0.172 

0.35 

0.24 

0.10 

1.67 

0.63 

S 

11 

J212928.92+000415.0 

0.180 

0.38 

0.30 

0.09 

0.35 

0.84 

E 

12 

J133323.14+503028.2 

0.197 

0.52 

0.51 

0.01 

0.31 

0.37 

S 

13 

J073313.08+313954.6 

0.198 

0.38 

0.29 

0.09 

1.15 

0.53 

S0/E 

14 

J073133.58+314113.0 

0.210 

0.48 

0.44 

0.05 

0.19 

0.73 

E 

15 

J172004.43+262701.2 

0.228 

1.05 

0.95 

0.10 

2.22 

1.6 

E 

16 

J080128.61+355046.1“’ 6 

0.231 

2.00 

1.84 

0.16 

2.82 

0.83 

I 

17 

J142511.92+374729.5 

0.233 

1.13 

1.07 

0.06 

1.06 

0.71 

S0/E 

18 

J084213.26+363020.5 

0.243 

1.17 

1.13 

0.04 

0.90 

0.54 

S 

19 

J131130.41-013216.1 

0.244 

0.99 

0.85 

0.14 

0.70 

0.64 

E 

20 

J084339.40+292025.2 

0.248 

0.97 

0.92 

0.05 

0.18 

0.34 

E 

21 

J082644.54+040705.4 

0.262 

1.59 

1.48 

0.11 

2.20 

1.2 

S 

22 

J 171957.79+264027.3 0 ’ 6 

0.263 

1.75 

0.59 

1.17 

8.11 

1.7 

E 

23 

J 133428.33+502829.0 

0.266 

0.62 

0.47 

0.15 

1.69 

1.1 


24 

J101641.15+384703.4 

0.269 

1.84 

1.78 

0.06 

1.13 

0.48 

S 

25 

J084209.68+293836.1 

0.279 

0.83 

0.74 

0.09 

1.91 

1.5 

S 

26 

J 133553.40+405459.2 

0.282 

1.32 

1.18 

0.15 

0.59 

1.5 

s 

27 

J 101653.99+390530.9 

0.292 

1.03 

1.00 

0.02 

0.43 

1.0 

s 

28 

J212914.75+001947.6 

0.306 

0.92 

0.86 

0.06 

1.09 

1.3 


29 

J073322.45+313915.5 b 

0.307 

0.99 

0.80 

0.19 

0.99 

1.1 

E 

30 

J010625.81+005343.3 

0.313 

1.09 

1.03 

0.06 

1.30 

1.6 

E 

31 

J010658.45+010146.8 

0.314 

0.63 

0.52 

0.11 

0.44 

0.45 

E 

32 

J010658.95+010438.3°’ 6 ’ c 

0.327 

0.95 

0.82 

0.13 

2.35 

0.66 

S0/E 

33 

J090016.83+205502.9 

0.333 

3.27 

3.06 

0.21 

2.44 

1.6 

E 

34 

J090034.67+204013.2°’ 6 ’ c 

0.352 

1.59 

1.02 

0.57 

12.53 

1.3 

E 

35 

J101805.93+385755.8°’ c 

0.369 

10.15 

9.64 

0.51 

8.89 


E 

36 

J213007.49+001419.1 a ’ 6,c 

0.395 

2.37 

2.02 

0.35 

4.12 

1.2 

S0/E 

37 

J 101742.86+385540.9*’ 

0.407 

1.32 

1.04 

0.28 

1.13 

0.84 

E 

38 

J024858.24-032446.9 a,i>,c 

0.429 

2.45 

1.86 

0.59 

2.43 

1.0 


39 

J133241.05+502502.9 

0.440 

4.76 

4.23 

0.53 

2.15 

3.2 

E 

40 

J101800.18+385833.5“’ 6,c 

0.440 

5.84 

3.83 

2.00 

8.22 

1.4 

I 
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Table 8—Continued 


# 

Source 

LIRAS 

z 

-^tot,IR 

1O 11 L 0 

Lsb.ir 

10" l. 0 

Lagn.ir 

1O 11 L 0 

^AGN, total 

1O 11 L 0 

0 

o 

Type d 

41 

J083759.22+145557.1 

0.456 

2.57 

2.25 

0.32 

1.31 

1.0 

E 

42 

J083244.18+654251.5°’ 6,c 

0.457 

4.33 

1.97 

2.35 

9.60 

0.92 

C 

43 

J073258.91+313724.5 

0.482 

3.65 

3.54 

0.11 

1.98 

0.62 


44 

J171919.29+262835.3 

0.507 

9.70 

9.27 

0.42 

2.87 

1.9 

E 

45 

J080049.66+360514.2°’ 6 ’ c 

0.511 

3.95 

3.26 

0.69 

12.12 

1.0 

C 

46 

J 133616.50+405529.4 

0.530 

2.38 

1.51 

0.87 

3.11 

0.60 

c 

47 

J101714.14+390124.4°’ 6,c 

0.536 

1.96 

0.89 

1.07 

18.72 

2.8 

I 

48 

J101600.54+391049.3 

0.538 

4.88 

4.78 

0.10 

2.07 

2.6 

I 

49 

J213011.84+000558.3°’ 6,c 

0.561 

4.62 

2.33 

2.28 

9.35 

1.5 

E 

50 

J082927.84+654906.5“’ 6,c 

0.568 

4.14 

3.98 

0.16 

2.76 

0.74 

C 

51 

J145635.24+222400. 9 a ’ b ’° 

0.590 

3.95 

3.09 

0.86 

4.77 

1.8 

c 

52 

J213015.48+000430.0“’ 6 ’ c 

0.604 

7.43 

6.57 

0.85 

4.43 

1.2 

c 

53 

J084317.56+293818.6°’ 6,c 

0.623 

5.21 

2.77 

2.44 

12.56 

1.3 

E 

54 

J015214.76+010705.7°’ 6,c 

0.702 

15.12 

12.96 

2.16 

11.22 

2.1 

c 

55 

J090126.15+205632.1 

0.756 

10.37 

10.26 

0.11 

1.95 

1.8 

E 


a High Luminosity Subsample (HLS) 
b AGN-dominated at 24 pm 
c Comparison sample (see text) 

d S = spiral, E = elliptical, C = too compact to classify, I = interacting 
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Flux(24 fji m) [mjy] 


Fig. 1.— 24 /im flux density distributions of sources in our survey area. The black histogram shows 
the distribution of all sources with 24 /.mi flux density > 1 mjy excluding stars and asteroids. The 
magenta curve shows the distribution of sources targeted by Hectospec fibers (1729). The green 
curve shows the distribution of sources with well determined redshifts from Hectospec spectroscopy 
(1209). The red histogram is Type-1 AGNs with emission line FWHM > 1200 km s -1 (205). 
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Fig. 2.— 24 /jrn flux density versus redshift for Type-1 AGNs with 24 flux density > 1 mjy in 
our survey area. The filled circles show Herschel -detected Type-1 AGNs. The unfilled circles show 
Herschel- non-detected Type-1 AGNs. 
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24/i.m flux density (mJy) 




Fig. 3.— (a): Simulation results showing the expected apparent K-band magnitude as a function of 
24 /irn flux for Type-1 AGNs, as functions of redshift (z = 0-3.6) and dust extinction ( Ay = 0-1.5). 
(b): Same as (a) but for r' band, (c): K-band magnitude (Vega) distribution of 205 Type-1 AGNs 
with 24 fim flux density > 1 mJy in our survey area, (d): Distribution of r' band magnitudes for 
targets that were put on Hectospec fibers and targets that are successfully identified with emission 
lines. 
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10 “’ 10 ° 10 ' 10 2 10 3 
Observed Wavelength (/xm) 



1CT’ 10° 10’ I0 2 10 3 

Observed Wavelength (/xm) 



10 - ’ 10 ° 10 ' 10 2 10 3 
Observed Wavelength (/xm) 


Fig. 4.— Examples of SEDs and decomposition results. The diamond points are the average fluxes 
at FUV, NUV, and five SDSS bands (u 7 , g 7 , r 7 , i 7 , z 7 ), J, K, and three WISE bands (3.4 ^xrn, 4.6 
/xm, and 12 /xm), 24 /jm, and five Herschel bands (100 /xm, 160 /xm, 250 /xm, 350 /xm, and 500 
/xrn). The solid lines show the SED decomposition results: the magenta line is the rescaled Type-1 
AGN template (Elvis et al. 1994); the green line is the stellar photospheric component; the blue 
line is the best fitted star formation template. The black line is the total of AGN, stellar, and star 
formation components. 
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10 “' 10 ° 10 ' 10 2 10 3 
Observed wavelength (fj, m) 



Observed wavelength Quin) 



Observed wavelength {fi m) 


Fig. 5.— The average SEDs for Type-1 AGNs with no formal FIR detections in our sample in three 
discrete redshift bins: z = 0.1 — 0.7 (top, 24 sources), 0.7 — 1.2 (middle, 24 sources), and 1.2 — 1.9 
(bottom, 18 sources). Symbols are the same as in Figure |4j 












-64- 



0.0 0.5 1.0 1.5 2.0 2.5 5.0 

Redshift 


Fig. 6.— Observed-frame [250/24 /mi] flux ratio versus source redshift for the Type-1 AGNs in 
our sample. The downward pointing arrows indicate the sources not detected at 250 jun, based on 
3a upper limits. The five-pointed stars are the average values for Herschel non-detected Type-1 


AGNs from our stacking analyses (see Section 

3.2.2 

). The dotted (red) line is the [250/24 /im] flux 

ratio of Type-1 AGN template from Elvis et al. 

1994) 

. The dashed (blue) line is the flux ratio of 

a typical Type-1 AGN template from 

Fritz et al. 

to 

o 

o 

rft 

)• 
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Fig. 7.— The FIR dust temperature versus total infrared luminosity (rest-frame 8 yum-1000 /im). 
The big circles (blue) show the L-T relation of star-forming galaxies (z=0.1-2) derived from Her- 
MES and PEP data in COSMOS, GOODS-S and GOODS-N fields (Symeonidis et ah, 2013). The 
error bars show the 1-sigma scatter of the L-T relation. The blue irregular polygon encloses the 
sources with strong warm infrared components. The remaining sources are compatible with the 
L-T relation, particularly if one allows for modest warm components in a few of them. 
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IQ-' 10 ° 10’ 10 2 10 3 

Observed Wavelength (/x m) 



ICT’ 10° 10' 10 2 10 3 

Observed Wavelength (/xm) 


Fig. 8.— SED decomposition results for an AGN at z = 1.883. Top: Decomposition with stellar, 
AGN, and star forming galaxy templates only. The measurements from 10 /xm to 200 ^m (3 /xm to 
60 /xm in the rest-frame) are high, indicating a warm excess from the MIR to FIR. Bottom: The 
same as the figure on the top. An additional warm component based on a model of a circumnuclear 
starburst has been added to improve the fit. 
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A (g'-i') 


Fig. 9.— (a): The observed-frame color ( g' — i') of Type-1 AGNs in our sample as a function of 
redshift. The filled circles (red) show Herschel -detected Type-1 AGNs; the unfilled circles (blue) 
show Herschel- non-detected Type-1 AGNs. The small dots (grey) represent the SDSS optically 
selected Type-1 quasars from the SDSS Data Release 7 Quasar Catalog (Schneider et al. 2010). 
The solid black line is the median value of the color of SDSS optically selected Type-1 quasars, 
(b): The K-S test shows that the relative color distributions of Herschel -detected (red solid) and 
-non-detected (blue dashed) AGNs are not statistically distinguishable (P-value=0.973). 
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Fig. 10.— Comparison of Type-1 AGNs in our sample and SDSS optically-selected Type-1 AGNs. 
(a): Redshift distribution. The red-solid histogram is Type-1 AGNs in our sample, while the black- 
dashed histogram is the SDSS optically-selected sample. The blue hatched histogram is the SDSS 
AGNs that are not included in our sample due to their 24 /mi flux density being below 1 mJy. (b): 
i' magnitude distribution. The symbols are the same as in the top panel. 
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Fig. 11.— Comparison of Type-1 AGNs in our sample and SDSS optically-selected Type-1 AGNs. 
(a): The relative observed-frame color A (g' — i') as a function of redshift. The red filled circles 
show Type-1 AGNs in our sample. The blue circles show the SDSS optically-selected AGN sample 
within the same area; they are open if the source is not in our sample. The small dots (grey) 
represent the SDSS optically selected Type-1 quasars from SDSS Data Release 7 Quasar Catalog 
(jSchneider et al. 2010). (b): The [24 /irn/2] flux ratio as a function of redshift. Symbols are the 
same as in the upper panel. The dotted and dashed lines are the flux ratio calculated from Elvis’s 
quasar template (1994) with reddening Ay = 0 and 0.5, respectively. 
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Fig. 12.— Upper: Emission-line diagnostic diagram for sources at z < 0.34 taken from Kewley et 


al. (2001} (Equation [ 2 ]) . Filled circles are Herschel -detected Type 2 AGNs. Unfilled circles are 


Herschel- non-detected Type 2 AGNs. Lower: The distribution of [O ill]A5007/H/3 as a function of 
redshift for AGNs with z > 0.34. 
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Observed Wavelength (/x m) 



Observed Wavelength (/xm) 


Fig. 13.— Examples of SED decomposition fits. Upper: SED dominated at 24 fim by the AGN. 
Lower: SED dominated at 24 fim by star formation. The cyan, green, magenta, and blue solid 
lines represent the best-fitting young stellar component, old stellar component, AGN component, 
and starburst component, respectively. The black solid line represents the total of the best-fitting 
models. 
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Observed Wavelength (/u.m) 



Fig. 14.— The average SEDs for Type-2 AGNs with no formal FIR detections in our sample in 
two discrete redshift bins: z =0.0-0.4 (upper, 14 sources), and 0.4-0.8 (lower, 13 sources). Symbols 
and line colors are the same as in Figure [13} 
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L[OIII] (*-©) 


Fig. 15.— The relation between AGN total luminosity and [O hi] A5007 luminosity for Herschel- 
detected Type-2 AGNs in our sample. The fitted line (all points with equal weights) has a slope of 
0.74. 
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Redshift 


Fig. 16.- The ratio of AGN luminosity to Eddington luminosity for our Type-2 AGNs. Filled 
circles are for AGN-dominated sources, while open ones are for SF-dominated. 
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(M„/M 0 ) (SDSS KCORRECT) 


Fig. 17.— Comparison of the stellar mass calculated from Bell et al. (2003) and SDSS KCORRECT 
(private communication, Krystal Tyler). The SDSS KCORRECT stellar masses are based on the 
Bruzual-Charlot stellar evolution synthesis. The stellar masses derived from Bell et al. (2003) are 
consistent with those from SDSS KCORRECT. 
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Fig. 18.— Simulation results for the ability of SED decomposition to constrain the stellar 
component in the NIR. The simulation is performed 10000 times for each fixed value of AGN 
and stellar component in the rest-frame J-band. The derived results for the input flux ratio 
fluxAGN,j/huxgteiiar,j = 0.5, 1, 2, and 3, are shown in red dash-dot, black solid, magenta dotted, 
and blue dashed lines. 
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Fig. 19.— Comparison of the stellar masses estimated by K-band luminosity using the equation 
from Bell et al. (2003) (See Section 6.4) and the stellar masses derived from the local mass ratio 
M*/M. = 700. The filled and unfilled circles are AGNs at z < 0.6 and 0.6 < z < 1.2, respectively. 
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A. Images of Type-2 Host Galaxies 
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Fig. 20.— Subaru images of Herschel -detected Type-2 AGNs. The circle radius is 15 w . The 1st 
part of a continued figure. 
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J40 J*' ^ #50 #51 #52 




Fig. 20.— Subaru images of Herschel -detected Type-2 AGNs. The circle radius is 15". The 2nd 
part of a continued figure. 


B. 24 /xm-Selected Type-1 AGN Sample in the LoCuSS Fields 

In total, we detected 2439 sources with 24 /jin flux above 1 mJy in the LoCuSS fields. Out of 
these, the following 541 sources were not included in the target list for the Hectospec spectroscopic 
follow-up: 

1. 71 sources that were outside the available near-infrared images. 

2. 168 sources that were identified as stars. 

3. 373 sources with no obvious optical/near-infrared counterparts (The 5-cr detection limit of 
the Subaru images in r or i band is ~ 25 magnitude and the 5-cr detection limit at K-band is 
19 mag (Vega)).). 

The remaining 1827 24 ^rnr sources are likely to be extragalactic. We may have discarded a number 
of extragalactic sources with faint optical/near-infrared counterparts (category 3 above) although 
some fraction of the category 3 sources is expected to be asteroids. 

Among these 1827 sources, 1729 were observed by Hectospec while another 18 sources have 
spectroscopic information from SDSS. The completeness of the spectroscopic coverage is therefore 
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about 94.6%. Among the 1729 sources targeted by Hectospec, 1263 sources have produced spec¬ 
troscopic redshifts with the corresponding success rate of 73%. However, the sources that did not 
produce spectroscopic redshifts are unlikely to be Type-1 AGNs. Therefore, our 24 /um-selected 
Type-1 AGN sample is expected to be complete at the ~94% level, which is the completeness of our 
spectroscopic coverage. Thus, we have 205 sources that satisfy our Type-1 AGN selection criteria 
(See Section [3]), 177 confirmed with Hectospec spectra, and 28 confirmed with SDSS spectra. 


C. Black Hole Mass Estimate 

The following methods were used to estimate black hole masses from our spectra: 


1. FWHM(H/3) and L a (5100A). For the optical continuum luminosity and FWHM of the 
H/3 broad component, 


log M B h(H/3) = log 


/FWHM(H/3)\ 2 /AT a ( 5100A)\ ( 
\1000 km s -1 / \ 10 44 erg s _1 J 


+ (6.91 ±0.02). (Cl) 


The sample standard deviation of the weighted average zeropoint offset is ±0.43 dex (Vester¬ 
gaard Sz Peterson 2006). 


2. FWHM(Mgll). For a given wavelength, A, the black hole mass based on Mg II was obtained 
according to: 

^ ^ (C2) 


where zp{\) is 6.72, 6.79, 6.86, and 6.96 for A1350A, A2100A, A3000A, and A5100A, 


re¬ 


spectively. The 1 a scatter in the absolute zero-points, zp , is 0.55 dex (Vestergaard &; Osmer 
2009). 


3. FWHM(CIV) and L a (1350A). For the ultraviolet continuum luminosity and the FWHM 
of the C IV line, 


log Mbh(Civ) = log 


/FWHM(Civ)\ 2 / AL a (1350 A) A 
\ 1000 km s _1 ) \ 10 44 erg s _1 J 


0 . 53 ' 


± (6.66 ±0.01). (C3) 


The sample standard deviation of the weighted average zeropoint offset is ±0.36 dex (Vester¬ 


gaard & Peterson 2006). The L A (1450 A) luminosity is equivalent to L A (1350 A) in the equa¬ 


tion above without error or penalty in precision (Vestergaard fe Peterson][2006). 


Figure 21 shows the examples of broad emission line fits for H/3, Mgn, and Civ, respectively. 


































- 81 - 



6000 6200 6400 6600 6800 
Observed Wavelength (A) 




4250 4750 5250 5750 6250 

Observed Wavelength (A) 


Fig. 21.— Examples of broad emission line fits, (a): H/3; (b): Mgn; (c): Civ. For each panel, 
the upper black line shows the original SED. The lower blue line shows the continuum and Fe 
subtracted SED. The upper magenta line shows the full fits; the lower magenta line shows the fits 
for the emission lines; the gray lines show the flux density errors; the green lines show the broad 
Gaussian components, while the red lines show the narrow Gaussian components. 


D. Correction of the AGN Template for Star Formation 


The template we use for the intrinsic AGN SED was built from a detailed set of observations 


of a representative set of optically selected quasars by Elvis et al. (1994). A more recent study 


by Richards et al. (2006) used a similar approach and derived a virtually identical template. The 


excellent agreement is encouraging; for example, our results are independent of which template 
we use. However, neither study attempted to correct the templates for the far infrared emission 
due to star formation. Doing so is challenging because one needs an independent, extinction-free 
estimate of the rate of star formation in the quasar host galaxies. The 11.3 /rrn aromatic feature 
is an appropriate indicator, particularly since it is not strongly affected by an AGN (Diamond 
Stanic & Rieke 2010). We have therefore used a large set of measurements of this feature in quasar 
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spectra (Shi et al. 2014), along with a star forming galaxy far infrared template (Rieke et al. 2009) 
to estimate the necessary correction. The approach was to correlate the equivalent width of the 
11.3 fim feature with the ratio of fluxes at 25 and 60 /trn (IRAS) or at 24 and 70 /trn (MIPS) to 
determine the influence of star formation on the far infrared spectrum in a variety of galaxies with 
and without AGN. We then used the relation derived from this correlation analysis and the average 


EW of the 11.3 /jm feature for the quasar sample used by Elvis et al. (1994) to determine how to 
adjust their template in the far infrared. 


The initial template we used was for radio-quiet quasars; Elvis et al. (1994) list 19 of these 


sources with IRAS detections, and they would have been most influential in determining the far 
infrared behavior of their template (we return to the IRAS upper limits later). Of those 19, we have 
11.3 n m EW measurements for 15 (79%), with an average value of 0.037 /im (standard deviation 
of the mean = 0.007 gm). A linear fit to the dependence of EWs vs. infrared flux ratios indicates 
that the ratio of IRAS 60 to 25 /iin flux densities for the Elvis et al. (1994) template has been 
boosted by a factor of 1.24 due to star formation, relative to the case for an EW of 0.0. However, 


the baseline in EW is small, so we repeated the determination adding the galaxies from Brandi 


et al. (2006) (which we selected because the methodology for determining EWs was similar to 
the method for the quasars). This reference includes cases with EW up to ~ 0.9, thus extending 
the baseline and improving the determination of the slope of the relation. This fit indicated a 
star-formation induced boost in the far infrared flux ratio for the Elvis template by a factor of 1.27. 
When we added the radio loud quasars in the Elvis sample plus additional PG quasars with 11.3 
/iin and far infrared measurements, and substituted MIPS for IRAS measurements when they were 
available, we got a value of 1.27. This last correlation is illustrated in Figure 21. 

With a determination of the size of the star-formation boost in the flux ratio, we subtracted 
a star-forming galaxy template (specifically for L(TIR) = 10 11 Lq (Rieke et al. 2009)) from the 


Elvis AGN template. We used synthetic photometry on the f60/f25 flux density ratio to match the 
results from the correlation analysis based on the EW of the 11.3 /mi feature. 

The adjusted AGN template may be an extreme case, since we did not include the galaxies 


in the sample of Elvis et al. (1994) for which there were only IRAS upper limits. These galaxies 


should include those with the weakest star formation relative to the AGN, as well as some that are 
just fainter then the detected ones at all wavelengths. It is not possible to reconstruct exactly what 
effect the upper limit cases would have had on the published template, but presumably they tended 
to make it fainter in the far infrared than it would have been based only on the IRAS detected cases. 
Thus, we consider our adjusted AGN template to be a limiting case for the maximum plausible 
far infrared contribution from star formation, and take the unadjusted template to be the limiting 
case in the other direction. 

This approach provides a correction out to 100 gm (rest). Beyond this wavelength, the Elvis 
template is a power law interpolation to the radio regime. There are very few examples of quasars 
that can be shown to have very low leves of star formation and at the same time have sufficiently 
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EW 11.3fj.m 


Fig. 22.— The relation between the equivalent width of the 11.3 /im aromatic feature and the 
ratio of flux densities at 70 and 24 fim (from MIPS) or at 60 and 25 ^m (from IRAS, if MIPS 


measurements are not available). The data are from Shi et al. (2014) and Brandi et al. (2006) 


sensitive measurements of upper limits at wavelengths longer than 100 /zm. Two examples, PG 
1501+106 and PG 1411+442, indicate that the power law substantially overestimates the fluxes 
in this region. Therefore, a more realistic replacement is a blackbody of 118K, with a wavelength 
dependent emissivity proportional to A -15 and scaled to match smoothly to the corrected SED at 
wavelengths short of 100 /zm. 












